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ABSTRACT 


This  study  is  concerned  with  the  derivation  of  approximate 
distribution  functions  for  sample  serial  correlation  coefficients  in 
linear  Markov  processes.  Two  cases  are  considered:  first,  when  the 
defined  process  is  circularly  correlated,  and  second,  when  the  process 
is  stationary  to  the  second  order.  Chapter  one  contains  a  review  of 
the  literature  relevant  to  this  special  problem.  In  Chapter  two, 
certain  results,  preliminary  to  the  application  of  asymptotic  procedures, 
are  established.  Asymptotic  expansions  for  the  distribution  functions 
are  derived  in  Appendix  three,  following  the  techniques  developed  by 
A.  Erdelyi  and  M.  Wyman,  and  in  Chapter  three  these  results  are  used 
to  obtain  asymptotic  representations  for  the  distribution  functions. 
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CHAPTER  I 

INTRODUCTION 


Even  a  cursory  survey  of  statistical  literature  suffices  to 
indicate  that  the  distribution  problem  for  serial  correlation  coefficients 
is  by  no  means  of  recent  origin.  A  considerable  volume  of  work  dealing 
with  various  aspects  of  this  problem  has  been  published  since  19^0,  and 
as  a  continuation  of  this,  we,  in  this  thesis,  deal  with  one  specialized 
aspect  of  this  problem.  Specifically,  we  concern  ourselves  with  the 
derivation,  by  asymptotic  procedures,  of  approximate  expressions  for 
distribution  functions  for  sample  estimates  of  serial  correlation  coeffic¬ 
ients  in  linear  Markov  processes. 

The  first  application  of  asymptotic  methods  to  this  problem  was 
made  in  DANIELS  (1956),  where  Debye's  method  of  steepest  descents  was 
used  to  obtain  approximations  for  probability  density  functions.  All 
prior,  and  a  majority  of  subsequent,  publications  deal  with  variations  of 
this  problem  by  non-asymptotic  procedures.  Since  in  the  ensuing  discussion 
we  shall  only  be  concerned  with  this  specialized  point  of  view,  that  is, 
the  application  of  asymptotic  procedures  to  the  distribution  problem,  we 
forego  a  survey  of  the  literature  in  the  non-asymptotic  setting,  referring 
the  interested  reader  to  PATTON  (I96I)  and  ABBASSI  (196k)  for  comprehensive 
summaries . 

We  now  give  a  brief  resume  of  the  development  in  DANIELS  (1956). 

Consider  an  ordered  set  (n-  tuple)  of  random  variables , which  is 
assumed  to  be  a  sample,  drawn  at  times  t^,...,tn,  from  the  linear  Markov 


2 


process  defined  by 


x 


s 


where  s  =  l,...,n,  a  is  a  real  constant,  called  the  lag  one  serial 


correlation  coefficient  with  |a|  <  1.  The 


are  independent 


random  variables,  each  having  a  normal  distribution  with  mean  zero  and 
variance  one.  There  is  imposed  the  further  restriction  that  the  state  of 
the  process  at  any  given  time  be  independent  of  assumptions  regarding  an 
initial  state.  The  devices  used  by  Daniels  to  realize  this  restriction 
are  two:  first,  following  a  suggestion  by  Hotelling,  he  assume  that  the 
process  defined  is  circularly  orientated,  that  is,  the  sample  observation 
at  time  t  1  is  taken  to  be  identical  with  the  sample  observation  at 


time  t^,  for  any  predetermined  finite  integer  n.  This  artificial 
restriction  not  only  obviates  any  necessity  for  the  consideration  of  a 
possible  initial  state,  but  also  simplifies  the  subsequent  analysis  of 
the  distribution  problem.  In  the  second  case,  Daniels  assumes  that  the 
process  defined  is  stationary  to  the  second  order  [which,  by  virtue  of 
the  normality  of  the  e's,  implies  complete  stationarity ] .  This 
assumption,  in  contrast  to  the  circularity  restriction,  is  a  natural  one, 
if  it  is  assumed  that  the  process  has  been  in  progress  for  a  considerable 
period  of  time. 


In  the  circular  case,  the  joint  distribution  of  the  x's  is 


multivariate  normal,  and  is  given  by 


dF(x  , ...,x ) 
i  n 


(2*)2 


dx. ...dx  , 

1  n 


and  the  maximum  likelihood  estimator 


3  - 


a 


x.x.+x.x  +..,+x  . x  +x  x. 

-I  2 — 2_J n-1  n n_L 


xl  + 


+  x£ 


is  used  as  the  sample  estimate  of  a  . 

In  the  non-circular  case,  that  is,  in  the  stationary  case,  the 
joint  distribution  of  the  x's  is  also  multivariate  normal,  and  is  given 
by 

dF(x1,...,xn) 

=  ■^■r-a“n  exp[-J{xJ+(l+a2)(x|+..„+x^_1)+x^-aa(x1x2+...+xn_1xn))] 

(at)5 

dx. ...dx  , 

1  n 


and,  the  intra-class  correlation  coefficients  between  the  sequences 
[x^ > • • • » xn_ ^ )  a°d  { x2 , . . , , x^ }  , 


a  = 


w 


. .  +x  ,x 
n-1  n 


12  2  2 

5Xf+X^+ ...  +  X.  . 
^12  n-1 


+^x2 

^  n 


is  used  as  sample  estimate  of  a  . 


By  a  consideration  of  the  joint  moment-generating  function  of 
Cq  and  C^,  together  with  an  application  of  the  Cramer— 3aary  inversion 
formula,  Daniels  derives  approximations  for  density  function  of  the  ratios 
a  =  Cl/C0  in  each  of  the  above  cases.  Explicitly,  the  expressions 
extracted  are,  for  the  circular  case 


-  k  - 


h(r)  [=  Pr(o=r)] 


r(^  +  ( 1-r2) 


n  _1 

2  "  2 


>/*  (§  +  |) 


( l-2ar+a2) 


and  for  the  stationary  case, 


h(r)  [=  Pr(o=r ) ] 


F(g  -1)  (l-q^2  (1-r^2 


5  -  1 


2  /*  r(5  -  i) 


( l-ar)(  l-2ar+a2)‘ 


?  -  1 


where,  in  both  cases,  -1  <  r  <  +  1, 


The  two  estimators  considered  are  ratios  of  quadratic  forms  in 
with  the  numerator  indefinite  and  the  denominator  positive 
definite,  which  suggests  the  consideration  of  the  more  general  problem  of 
obtaining  expressions  (approximate  or  otherwise)  for  distributions  of 
ratios  of  quadratic  forms,  with  indefinite  numerators  and  positive  definite 
denominators . 


In  GURLAND  (19^8),  ( 1953) >  (1955)  and  (1956)  are  found  various 
aspects  of  this  problem  in  its  more  general  setting.  The  approach  taken 
in  these  papers  is  briefly  as  follows:  and  are  taken  to  be 

arbitrary  quadratic  forms  in  the  random  variables  x^,...,xn.  The 
x, ,...,x  have  a  joint  multivariate  normal  distribution,  and  is 

restricted  so  that  it  is  always  positive  definite.  The  linear  function 
u  =  -  tCq  (linear  homogeneous  in  C^,  with  -1  <  r  <  l)  is 

considered  next,  which  by  the  manner  of  its  definition  is  an  indefinite 
quadratic  form  in  x^,...,xn.  Then,  in  contrast  to  DANIELS  (1956),  the 
characteristic  function  of  u  is  considered,  and  an  inversion  formula 
for  this  characteristic  function  is  derived  [cf  GURLAND  (I9I4.8)  and  GIL- 
PALAEZ  (I95I)].  The  distribution  function  of  the  ratio  C^/C^  is  then  derived 


. 


■ •••  '  •  .  ■  .  .  . : ,  -  '  •  .  ■■ : 
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by  setting  u  =  0  in  this  inversion  formula.  Thus,  effectively,  the 
problem  in  its  most  fundamental  form  is  that  of  the  derivation  of  distri¬ 
butions  of  indefinite  quadratic  forms  in  x, ,...,x  . 

1  n 

In  GURLAND  (1953)  and  (1955)»  the  procedure  outlined  above  is 
utilized,  and,  convergent  Laguerrian  expansions  for  distribution  functions 
for  quadratic  forms  [positive  definite  and  indefinite],  and  ratios  of 
quadratic  forms  [denominator  always  positive  definite]  are  derived.  In 
GURLAND  (1956)  the  same  results  are  modified,  in  that  convergent  expansions 
for  quadratic  forms  and  ratios  of  quadratic  forms  [with  the  same  restriction 
as  before]  are  derived  in  terms  of  Incomplete  Gamma  and  Incomplete  Beta 
functions  respectively. 

In  this  work,  our  approach  is  almost  exactly  that  of  DANIELS 
(1956).  That  is,  we  assume  the  existence  of  linear  Markov  processes  under 
the  restrictions  of  circularity  and  stationarity  assumptions,  and,  further, 
use  as  estimators  of  the  serial  correlation  coefficients,  the  maximum 
likelihood  estimator  and  the  intra-class  serial  correlation  coefficient 
respectively.  However,  in  contrast  to  DANIELS  (1956),  we  derive  approx¬ 
imations  for  distribution  function  rather  than  density  functions. 

Notwithstanding  the  restricted  nature  of  our  work  in  comparison 
with  that  of  Gurland,  namely,  that  of  considering  special  sample  serial 
correlation  coefficients  rather  than  more  general  ratios  of  quadratic 
forms,  and  the  difference  in  our  terminal  objectives,  namely,  that  of 
obtaining  asymptotic  expressions  rather  than  convergent  expansions,  our 
work  displays  definite  evidence  of  Gurland' s  influence,  especially  in 


Chapter  two. 


J  '  '.  0. 


' 

!/  • V.  >hci'i,  b  ■  •  .  ...  if.:-; 
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Briefly,  the  procedure  in  the  ensuing  chapters  and  appendices 
takes  the  following  direction. 

In  Chapter  two,  we  obtain  the  complex  characteristic  function 
for  u  =  C ^  ~  tCq  [where  C^,  take  the  meanings  assigned  with  respect 

to  sample  serial  correlation  coefficients  and,  as  before,  -1  <  r  <  +  l]. 
Following  techniques  similar  to  DAN IEI5  ( 1956) ,  we  utilize  the  analyticity 
properties  of  the  characteristic  function  [cf  WIDDER  (I9UI)  and  LUKACS 
(i960)]  to  derive  complex  integrals  which  represent  the  distribution 
function  for  the  sample  serial  correlation  coefficient. 

In  Appendix  three,  diverging  from  Daniels'  approach,  we  use 
techniques  developed  in  ERDELYI  and  WYMAN  (I963)  and  (1963a)  to  derive 
asymptotic  expansions  for  the  complex  integral  in  terms  of  Confluent 
Hypergeometric  functions. 

Finally,  in  Chapter  three,  using  the  results  obtained  in 
Appendix  three,  we  give  Asymptotic  representations  for  distribution 
functions  of  the  defined  sample  serial  correlation  coefficients. 

Appendices  one  and  two  are  devoted  to  certain  auxiliary  results 
used  in  the  developments  of  Chapter  two,  Appendix  three  and  Chapter  three. 

Appendix  four  contains  graphical  illustration  of  the  behavior 
of  the  derived  distribution  functions  for  various  values  of  a  and 


various  sample  sizes  n 


. 
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CHAPTER  II 

PRELIMINARY  RESULTS 

We  wish  to  derive,  by  Asymptotic  Procedures,  approximations  to 
distribution  functions  for  sample  estimates  of  lag  one,  serial  correlation 
coefficients.  The  stochastic  process  is,  specifically,  the  linear  Markov 
process  in  discrete  time  defined  by 

Xt  -  aXt-l  =  €t 

for  all  t  ,  where  the  constant  a  is  called  the  serial  correlation 
coefficient  and  the  (€t)  are  independently  and  identically  distributed, 
N(0,l)  ,  random  variables. 

§  1.  Two  Special  Cases: 

In  this  thesis  we  confine  our  attention  to  two  special  cases.  In 
the  first  case  we  assume  that  the  process  defined  is  circular,  that  is, 

(2.1.1)  xt-axt-l  *  £t 

for  t  =  1, . . . ,n  ,  where,  |a|  <  1  and  x^  =  .  In  the  second  case 

we  assume  the  process  to  be  stationary,  that  is, 

(2.1.2)  xfc  -  ccxt_1  =  €t 

for  t  =  l,...,n  ,  where  |a|  <  1  and 

/  \  cJS  I 

cov  (V  xt-s>  *  • 

1-a 


■ 


. 


-  -  -  c  •  ' 
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Let  x, be  a  sample  of  observations  from  either  one 
1'  1  n  r 

of  the  two  defined  processes.  They  have  a  Joint  multivariate  normal 
distribution  since  are  independent,  N(0,l),  random  variables. 


Explicitly,  for  the  circular  case,  the  distribution  of  x^, 
is  given  by 


•  t 


x 

n 


(2.1.5)  dF(x1,...,xn) 


(i-qn) 

(2*) 


^72  exp[-i((1-+q2)(xi+- 


2 

,  +xn)-2a(x1x2+x2x5+. .  •+xn_1xn+xnx1 ) )  ] 


dx.  . .  .dx  ^ 
1  n  ’ 


which,  using  the  left-hand  matrix  notation,  may  be  expressed  as 


(2.1.4)  dF(xx, . . .,xn)  =  exP[~|  x>  x]  dx1...dxn  , 


(2*)' 


where, 


£  -  (x^> • • • ) xn) 


and 


-  9  - 


1-fa2 

-a 

0 

0 

•  •  • 

0 

0 

0 

-a 

-a 

1-fa2 

-a 

0 

000 

0 

0 

0 

0 

0 

-a 

1-fa2 

-a 

000 

0 

0 

0 

0 

0 

0 

-a 

1-fa2 

000 

0 

0 

0 

0 

0 

0 

0 

0 

0  0  0 

1-fa2 

-a 

0 

0 

0 

0 

0 

0 

0  0  0 

-a 

1-fa2 

-a 

0 

0 

0 

0 

0 

0  0  0 

0 

-a 

1+a2 

-a 

-a 

0 

0 

0 

0  0  0 

0 

0 

-a 

1+a2 

with 


12" 1 F 


=  1  -  a 


n 


For  the  non-circular  case,  the  distribution  of  x, ,  ...,x  is 

’  V  7  n 


given  by 


(2.1.5)  dF(x1,...,xn) 


P  ! 

d-q  )2 

(2*) 


't/I  exp[-i(x1+(l«.  )(x2+...+xn_1)+xn-2a(x1x2+...+xn_1xn))] 


dx  . . .dx 
1  n 


which  expressed  in  left-hand  matrix  notation  is 


(2.1.6) 


2  -1 

dF(x  ,...,x  )  >  '1~a3p  exP[*i  x'  2"1  x]  dx-i 
1  °  (2x)n/2 


.dx 

n 


where, 


10 


and 


t 


1 


with 


1 

-a 

0 

0 

fit 

0 

0 

0 

0 

-a 

1-kl2 

-a 

0 

•  •  • 

0 

0 

0 

0 

0 

-a 

1-KL2 

-a 

•  •  • 

0 

0 

0 

0 

0 

0 

-a 

1-KL2 

•  •  • 

0 

0 

0 

0 

0 

0 

0 

0 

•  •  • 

1  2 
1-KL 

-a 

0 

0 

0 

0 

0 

0 

•  •  • 

-a 

1+a2 

-a 

0 

0 

0 

0 

0 

•  »  » 

0 

-a 

1-kl2 

-a 

0 

0 

0 

0 

•  •  • 

0 

0 

-a 

1 

IM"1^  =  (l-a2)2  . 


As  an  estimate  of  a  we  use,  in  the  circular  case,  the  function 


(2.1.7) 


A  / 

a(x1; 


n 


)  = 


x.x  +x_x  +...+X  .x  +x  x, 
12  2^  n-1  n  n  1 

2  2  0  2 
x,+x0+...+x  .  2+x 

1  2  n-1  n 


which  in  matrix  notation  is 


x'  N  x 


(2.1.8) 


A 

a 


x'  D  x 


where 


and 


N  = 


D  = 


0 

1 

"5 

0 

0 

•  •  • 

0 

0 

0 

1 

? 

1 

0 

1 

2 

0 

•  •  • 

0 

0 

0 

0 

0 

1 

0 

1 

2 

•  •  • 

0 

0 

0 

0 

0 

0 

1 

2 

0 

•  •  • 

0 

0 

0 

0 

0 

0 

0 

0 

0  0  0 

4 

0 

1 

2 

0 

0 

0 

0 

0 

0 

0  •  • 

1 

2 

0 

1 

2 

0 

0 

0 

0 

0 

9  0  0 

0 

1 

2 

0 

1 

2 

1 

2 

0 

0 

0 

0*4 

0 

0 

1 

2 

0 

1 

0 

0 

0 

0  0  0 

0 

0 

0 

0 

0 

1 

0 

0 

0  0  0 

0 

0 

0 

0 

0 

0 

1 

0 

0  0  9 

0 

0 

0 

0 

0 

0 

0 

1 

4  9  0 

0 

0 

0 

0 

0 

0 

0 

0 

9  0  0 

1 

0 

0 

0 

0 

0 

0 

0 

4  0  0 

0 

1 

0 

0 

0 

0 

0 

0 

•  «  1 

0 

0 

1 

0 

0 

0 

0 

0 

0  0  4 

0 

0 

0 

1 

and  in  the  non-circular  case,  the  function 


12 


(2.1.9) 


X1X2 


+  xr 


+ 


+ 


+  X 

n- 


1 


x 

n 


2 

+  x  .  + 

n-1 


which  in  matrix  notation  is 

N  x'  N  2S 

(2.1.10)  a  =  - 

x'  D  x 

where, 


0 

1 

2 

0 

0 

•  •  • 

0 

0 

0 

0 

1 

2 

0 

1 

2 

0 

0  0  0 

0 

0 

0 

0 

0 

1 

2 

0 

1 

2 

•  •  • 

0 

0 

0 

0 

0 

0 

1 

2 

0 

•  •  • 

0 

0 

0 

0 

0 

0 

0 

0 

•  »  t 

0 

1 

2 

0 

0 

0 

0 

0 

0 

•  •  • 

1 

2 

0 

1 

2 

0 

0 

0 

0 

0 

0  0* 

0 

1 

2 

0 

1 

2 

0 

0 

0 

0 

9  0  0 

0 

0 

1 

2 

0 

and, 


-  13  - 


— 

1 

£ 

0 

0 

0 

•  •  • 

0 

0 

0 

0 

0 

l 

0 

0 

•  •  • 

0 

0 

0 

0 

0 

0 

1 

0 

0  0  0 

0 

0 

0 

0 

0 

0 

0 

1 

0  0  0 

0 

0 

0 

0 

0 

0 

0 

0 

0  0  0 

1 

0 

0 

0 

0 

0 

0 

0 

0  0  0 

0 

1 

0 

0 

0 

0 

0 

0 

0  0  0 

0 

0 

1 

0 

0 

0 

0 

0 

0  0  0 

0 

0 

0 

1 

2 

In  both  cases 

|a|  <  1  . 

Thus,  the  problem  of  finding  the  sampling  distribution  of  a  is  effectively 
that  of  finding  the  distribution  function  for  the  ratio  of  two  quadratic 
forms . 

§  2.  The  Distribution  of  the  Ratio  of  two  Quadratic  forms. 

We  have  the  ratio 

x*  N  x 

(2.2.1)  a  =  — - 

x*  D  x 

where, 

x*  =  (x^,..., xn)  , 

N  is  an  indefinite  symmetric  matrix  of  constant  coefficients, 

D  is  a  positive  definite  symmetric  matrix  of  constant 
coefficients , 


. 


- 
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Denoting  the  distribution  function  of  a  by  !!(•)  ,  by  definition 


(2.2.2) 


H(r)  =  Pr[a  <  r] 


=  P 


r  x'  N  x 


-  x'  D  x 


<  r 


where,  r  is  any  real  number  satisfying  |r|  <  1  .  For  reasons  not 
immediately  obvious,  which  become  apparent  as  the  discussion  unfolds,  we 
restrict  r  so  that  |r|  <  1  .  Thus, 


Putting 


[x*  N  x 
x'  D  x 


<  r 


P  [(x»  N  x)  <  r(x'  D  x) ] 


=  P  [x'  (N  -  r  D)  x  <  0]  . 


x' (N  -  r  D)  x  =  u(x1,...,xn) 


and  denoting  the  distribution  function  of  u  by  G(»)  ,  we  get 

(2.2.3)  H(r)  =  Pr[a  <  r]  =  Pr[u  <  0]  =  G(0)  . 

(N  •  r  D)  is  a  symmetric  matrix  which  is  negative  definite,  singular  or 
positive  definite  according  as  the  value  of  r  varies  over  the  open 
interval  (-1,  l),  that  is,  (N  -  rD)  is  an  indefinite  matrix  of  variable 
coefficients . 


The 


characteristic  function  of  u(x, ,...,x  )  is 

\  1?  7  n/ 


00  r,00 


*„«)  ■  /  / 

j  -00  j  -00 


p.*MXn)]  •  •  *>Xn) 


(2.2.4) 


exp[  i^u(x, 


where,  £  is  a  complex  variable 


u(x_,...,x  )  =  x'(M  "  r]>)  2£ 


and 


cp^(^)  exists  for  all  £  such  that  Im  £  =  0  ,  but  may  also  exist  for 
some  £  such  that  Im  £  ^  0  .  In  the  former  case  Cp^(£)  is  merely  a 


complex  valued  function  of  a  real  variable  (Re  £)  defined  on  the  line 
Im  £  =  0  ,  but  in  the  latter  case  cp^(£)  i-s  an  analytic  function  which  exists 
and  is  regular  in  some  strip  of  non-zero  width,  which  contains  Im  £  =  0  . 
However,  since  r  in  x T (N  -  rD)  x  is  not  "integrated  out",  r  behaves 
as  a  "free  variable",  and  a  continuous  variation  of  r  in  (-1,1)  causes 
the  width  of  the  strip  to  fluctuate.  To  determine  the  distribution  function 
of  u  from  the  characteristic  function,  we  use  an  inversion  formula  which 
is  independent  of  this  fluctuation.  [Since  we  do  not  rule  out  the  possibility 
of  a  degenerate  strip  for  some  values  of  r].  The  inversion  formula  is 


which,  for  v  =  0  ,  will  give  the  distribution  function  of  a  ,  that  is 


(2.2.6) 


where, 


...  '■  -  '  .-y-i  ;  1 7.‘  s‘  .•  r,;Jbw  ?5it1 

■ 


(2.2.7) 


lire 

T  — ►  oo 

e  -+  o 
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) 


the  path  of  integration  being  the  line  Im  £  =  0  . 


In  matrix  notation  the  characteristic  function  is 


\IH  ^ 


r°°  r00  -i 

/  ...  /  exp[i(;  x'(N-rD)x]  expf-^x'Z  x]  dx.  ...dx 

J  -  00  J.  00  -  -  -  -  I  n 


121 


2  Ooo 


(2«) 


n/2 


...  r  exp[-|-x'{2  1-2i^(N-rD)  }x]  dx^...dxn  . 
-00  -00 


There  is  an  orthogonal  transformation  that  will  diagonalize  ^-2i£(N-rD)}  . 
Making  this  transformation  and  integrating  we  get 


(2.2.9)  q>u(5)  = 


12  I4 

- i - T 

|g"  -21£(N-rD) | 2 


where 

Ig"1  -  2i£(N-rD)|"^ 

is  the  Jacobian  of  the  transformation. 

§  Characteristic  Functions  for  Special  Cases. 

In  this  section  we  give  explicit  expressions  for  the  characteristic 
functions  in  the  two  special  cases  considered  in  this  thesis. 
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The  circular  case: 

We  have,  by  (2.1.4), 


(2.5.1)  g'1  = 


with 


and  because  of 


(2.3,2)  N  -rD 


1-KX2 

-a 

0 

•  •  • 

0 

0 

-a 

-a 

1-kx2 

-a 

•  •  • 

0 

0 

0 

0 

-a 

1-KX2 

•  00 

0 

0 

0 

0 

0 

0 

•  0  • 

1-kx2 

-a 

0 

0 

0 

0 

0  0  0 

-a 

1-KX2 

-a 

-a 

0 

0 

0  0  0 

0 

-a 

1-kx2 

— 

i  ■■  i 

Ig"1^  =  1  -  a” 

(2.1.8) 


-r 

1 

2 

0 

•  •  • 

0 

0 

1 

2 

1 

2 

~r 

1 

2 

•  •  • 

0 

0 

0 

0 

1 

2 

-r 

•  •  • 

0 

0 

0 

0 

0 

0 

•  •  • 

-r 

1 

2 

0 

0 

0 

0 

•  99 

1 

2 

-r 

1 

2 

1 

2 

0 

0 

•  •  0 

0 

1 

2 

-r 

--  — 

— 
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(2.3.3) 


=  (a+iOn 


Tints 


2"1  -  2i£  (N 

-  rD) 

l+a2+2ir£ 

-(a+iO 

0 

0 

0 

-(a+iO 

-(a+iO 

l-KX2+2ir£ 

-(a+iO 

0 

0 

0 

0 

-(a+iC) 

l+a2+2ir^  ... 

0 

0 

0 

0 

0 

0 

l+a2+2ir£ 

-(a+i£) 

0 

0 

0 

0 

-(a+iO 

l+a2+2ir£ 

"(a+i  0 

-(a+i^) 

0 

0 

0 

-(a+i?) 

l+a2+2ir^ 

l+a2+2ir^ 

-1 

0 

0 

0 

-1  | 

a+i£ 

-1 

l+a2+2ir£; 

a+i£ 

-1 

0 

0 

0 

0 

-1 

l*x2+2irC 

o  lc  9  9  9 

a+i£ 

0 

0 

0 

0 

•  •  • 

l+a2+2ir£ 

a+i£ 

-1 

0 

0 

-1 

l-Hx2+2ir^ 

a+i£ 

-1 

-1 

0 

-1 

l+a2+2irC 

a+if; 

(a  +  iOn  |B|  . 
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Thus,  by  equation  (2.2.9), 

(2.3.4) 

,  n  i 

„  / * \  1-a  1 

9u  ^  "  7  .  i.  \n/2  77T" 

(o+iO  /  |B|- 

Now  | B | 

is  a  circulent  and  has  value 

(2.3.5) 

/,  nv2 

|  B  |  =  ii^-L 

1  n 

z 

where, 


(2.3.6) 

1  l-KX2+2irC 

z  +  —  =  - rs - i 

z  a+i£ 

Therefore, 


(2.3.7) 

2 

, „  z(l-2ar+a  ) 

a  +  i^  - 

(l-2rz+z  ) 

and 


(2.3.8) 

„  \ \  l-an  (l-2rz+z2)n/2 

^P..  ( ^  ( z ) )  -  P  n/P  ^ 

(l-2ar-nx2)  /2  (l-zn) 

where. 


(2.3.8a) 

E(z)  =  ia  -  .i . 

l-2rz+z 

The  non-circular  case: 


We  have  by  (2.1.6) 
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(2.3.9)  £_1  = 


1 

-a 

0 

•  •  • 

0 

0 

0 

-a 

1-kl2 

-a 

•  •  • 

0 

0 

0 

0 

-a 

1-kx2 

•  •  • 

0 

0 

0 

0 

0 

0 

•  •  • 

1-kx2 

-a 

0 

0 

0 

0 

•  •  • 

-a 

1-KX2 

-a 

0 

0 

0 

•  •  # 

0 

-a 

1 

with 

I2'1 1  -  1-  ^ 

and  because  of  (2.1.10) 


(2.3.10)  N-rD= 


—— 

— 

_ir 

2l 

1 

2 

0 

...  0 

0 

0 

1 

2 

-r 

1 

2 

...  0 

0 

0 

0 

1 

2 

-r 

...  0 

0 

0 

0 

0 

0 

...  -r 

1 

2 

0 

0 

0 

0 

1. 

•  ♦  •  ^ 

-r 

1 

2 

0 

0 

0 

...  0 

1 

2 

_lr 

21 

Thus, 

£l  -  2i^(N-rD) 


(2.3.U) 
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(a+iOn 


1+16 

-(a+iO 

0 

0 

0 

0 

-(a+iC) 

l-Kx2+2ir^ 

■ (a+i£ ) 

0 

0 

0 

0 

-(a+iO 

l4a2+21rC 

...  0 

0 

0 

0 

0 

0 

...  l+a2+2ir(; 

“(a+16) 

0 

0 

0 

0 

T-l 

+ 

• 

• 

• 

l+a2+2ir£ 

“(a+i£) 

0 

0 

0 

0 

-(a+16) 

1+16 

a+i£ 

-1 

0 

...  0 

0 

0 

-1 

l-Kx2+2irt; 

a+i(; 

-1 

...  0 

0 

0 

0 

-1 

l+a2+2irC 

a+i£ 

...  0 

0 

0 

0 

0 

0 

l+a2+2irC 

•  •  •  .  if 

a+i£ 

-1 

0 

0 

0 

0 

-1 

l-Kx2+2irC 

a+it, 

-1 

0 

0 

0 

0 

-1 

1+H 

a+il«J 

(a+i£ )n  |B*| 
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Thus,  by  equation  (2.2.9) 


(2.3.12) 


<PU(C) 


(i-g2) 

(a+H)n/2 


Now  | B* |  is  a  continuant,  and  has  value 


(2.3.13)  |B*|  = 


( / i  \2 /  \2  2n-2  f  \2/  \2-> 

{(1-az;  ( l-ar+az-rz )  -  z _ (a-z)  (a-r+z-arz)  J 

zn(l-z2)(l-2ar-Kx2)2 


where. 


(2.3.14) 


Therefore, 


(2.3.15) 


1 

z  +  - 


z 


l+g2+2ir C 

a+it; 


(a+i  £ )  = 


z 


(l-2ar+a2) 

2 

l-2rz+z 


and 


(2.3.16)  q>  (£(z))  = 


p  JL 

(i-g  )a 


U 


(l-2ar+a2) 


--1 

2  l 


(l-z2)2  (l-2rz+z2 )n^2 


{(l-az)2(l-ar-Ktz-rz)2-z2n  2 (a-z)2 (a-r+z-arz)2}2 


where 


£(z)  =  ia  - 


iz(l-2ar+a  ) 

(l-2rz+z2) 


(2.3.16a) 
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Now, 


(2.3.17) 


c  (z) 

fttr 


(l-2ar4a^) (1-2^ ) 

p 

(l-2rz+z  ) (l-az) (z-a) 


Huis,  corresponding  to  (2.3.7)  we  have 


V  (C(z))5' (z) 

(2.3.18)  - ^ - 


1-a 


(l-2ar-Kx,2)2 


n 


(l-z2)  (l-2rz+z2)2 

(z-a) (l-az) (l-zn) 


-1 


and,  corresponding  to  (2.3.16)  we  have 


^..(C (z ) ) C  '  (z) 

(2.3.!9)  - ^ - 


(1-a2)2 


n 


(l-2ar-HX2)2 


-1 


(i  2\2 

. (i-£j _ 

(z-a) (l-az) 


-  -1 

2  2 

v  _ (l-2rz+z  ) _ _ 

( (l-az)2(l-ar+az-rz)2+z2n  2(a-z)2(a-r+z-arz)2  )2 


Since  we  are  prepared  to  tolerate  errors  of  0(An)  [ |A |  <  1]  in 

2 

the  final  approximation,  for  (1-r  )  not  small,  we  omit 


1-a 


1-z 


n 


in  (2.3.18) 


and  replace 


{ (l-az)2(l-ar+az-rz)2  + 


2n-2  /  \2  /  \  2 

z  (a-z)  (a-r+z-arz)  )d 


hy. 


(l-az) (l-ar+az-rz) 


in  (2.3.I9) 
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Thus,  neglecting  terms  which  are  exponentially  small,  (2.3*18)  and  (2.3.19) 
become  respectively 


(2.3*20) 


<pu(e(«m'(*) 


TFT 


—  -i 

?  ? 

(l-2ar+a  ) 


(l-22)(l-2rz+z2^2 

(z-a) (l-az. 


and 


(2.3.21) 


<pu(C(zK'(z) 

TF1 


p  1_ 

hi 


£  -2 

(l-2ar-nx2)2 


5  B  .! 

(l-z2)2(l-2rz+z2)2 

2 

(z-a)(l-az)  (l-ar+az-rz) 


§  4,  Discussion  of  a  Conformal  Transformation. 

In  the  last  section  we  saw  that  the  mapping 

l-Kx2+2ir^  _  1 

a+i£  ”  2  +  z 


l£P 

9U(0  =  - * - —  d^ 

I2"1  -  2ir5(H-rD)|? 


(2.U.1) 

trans formed 

(2.4.2) 


into  an  expression  of  the  form 

<p„ (e(OH'(*)  K,  2  s 

(S.k.3)  - -  dz  ~  Yz'-aJ  (1~2rz+z  )  dz 
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where. 


1)  £  =  £(z)  represents  the  transformation  (2.4.1), 

2)  K'  is  a  function  of  (a,r)  , 

5)  g(z)  is  an  algebraic  function  taking  real  values  for 
z  real  and  is  regular  in  a  simply-connected  subset 
of  |z|  <  1  ,  which  contains  z  =  a  and  z  =  r  , 

4)  s  is  a  real  variable  taking  large  values  [s  =  —  -1]  . 


[corresponding  to  (2.3.20) 


K*  =  (l-2ar-Kt2)  2 


and 


g(z)  =  (l-z2)/(l-az) 


and  to  (2.3.21) 


K' 


(l-2ar-Kx2)^ 


and 


3 


•] 


Equation  (2.4.1)  maps  the  £-plane  cut  from  -i<»  to 


and 
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from  0 r (  to  +  i«  onto  the  open  disc  |z|  <  1  ,  with  the  cut  being 

i  /  1  ^\2 

mapped  onto  the  unit-circle  |z|  =  1  .  The  point  £  =  -  r~J~  corresponds 

i  ^  1  *4cc  ^  * 

to  the  point  z  =  1  and  the  point  £  =  o7 i'+r y  corresponds  to  the  point 

z  =  -  1  .  Further,  £  =  -oo  +  It  and  £  =  +  oo  +  ir*  (t,  t*  arbitrary) 
are  mapped  onto  z 


=  r  -  hj  T-r^  and  z  =  r  +  hi l^r^  respectively. 


=  r  +  hfl-i  2 


The  z-image  of  Im  £  =  0  is  the  curve  given  by  the  equation 


-2  2--  2-  2  -2  2  - 

(2.4.4)  2<x(zz)  -(l+2ar-t<x  )zz(z+z)+4r (l+a  )zz+2a(z  +z  )-(l+2ar-fa  )(z+z)+2a  =  0  , 


which,  in  cartesian  coordinates  with  z  =  x  +  iy  ,  is 


(2.4.5)  a(x2+y2)2-(l+2ar+a^)x(x2+y2)+2r (l+a2) (x^+y2)+2a(x2-y2)-(l+2ar+a2)x+a  =  0 


This  curve,  which  we  shall  henceforth  denote  by  T  ,  is  symmetric  about  Im  z  =  0, 
which  is  normal  to  the  curve  at  the  point  of  intersection  z  =  a  [the  z-image 
of  the  point  £  =  0].  Thus,  points  below  Im  z  =  0  lying  on  T  are  complex 
conjugates  of  points  lying  on  the  curve  above  the  line  Im  z  =  0  .  T  runs 
from  the  point  z  =  r  -  ivl  [the  z-image  of  the  point  Re  £=  -  00]  to  the 
point  z  =  r  +  [the  z-image  of  the  point  Re  £  =  +00]. 

Thus  we  define  the  z-correspondent  of  the  ^-integral  in  the  inversion 
formula  (2.4.6)  to  be 


(2.4.6) 


g(z)  (l-2rz+z2)S  dz 


[ 


1 

2jri 


9U(?) 


d? 


where, 


-  27  - 


(2.4.7)  t 

J[r] 

The  integrals  defined  are  curvilinear  on  T  ,  a  and  its  complex  conjugate 
a  are  points  on  T  as  are  z*  and  z*  ,  and  the  notation 

E 
— > 

here,  and  in  the  subsequent  discussion  denotes  passage  to  the  limit  through 
the  set  E  .  Thus  the  z-integral  exists  in  the  sense  of  a  Cauchy  Principle 
value  as  defined  in  Appendix  I. 

Now  by  virtue  of  (2. 5. 7)  and  (2.3.16),  it  is  clear  that 
|  g  ( z )  |  <  co  for  z  =  r  +  is/l  -r2  . 


r  ra 

rz*  l 

=  lim 

i  * 

/ 

r 

-  z* 

-'a 

a-*  a 

z^  r+i*A-r2 

Thus 

(2.4.8)  lim  ~~  g(z)  (l-2rz+z2)S  =  0  , 

2  Z"a 

z-»  r+iv  1-r 


and 


(2.4.9) 


lim 


z  ~*r 


-  is/  l—i 


1 

z-a 


g(z)  (l-2rz+z2)S 


0 


Thus  we  write 


(2.4.10) 


lim 

r 

a  -» a 
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where  A 


r  + 


1 \/l  -r‘ 


§  5.  Deformation  of  the  Path  of  Integration, 

We  are  concerned  with  integrals  of  the  form 


(2.5.1) 


ifr  f[r]  T^ST  g(z)  (!-2rz«2)s  > 


where, 


-P  =  lira 

•'[r]  r 


r  +  T  1  ,  (A-rWC?) 

Jl  J  a  - 


a-+a 


and 


g(z)  is  regular  in  |z|  <  1  . 


The  integrand  in  (2.5.1)  has  critical  point  at  z  =  r  ,  so  that  we  deform 
the  path  of  integration  to  pass  through  z  =  r  .  Three  special  cases  have 
to  be  considered.  They  are,  r<a,  r=a,  r>a  .  For  the  present,  we 
shall  deal  with  r  <  a  and  r  >  a  only,  leaving  the  consideration  of 
r  =  a  to  be  dealt  with  at  a  later  stage. 


The  deformation  of  the  path  is  always  made  in  a  manner  such  that 
(2.5.2)  |1  -  2rz  +  z2j  <  (1  -  r2) 

for  all  z  on  the  deformed  path,  which  we  denote  by  L  .  Since  the  integrand 
is  regular  for  all  z  ,  except  z=a,in  |z|  <  1  ,  no  residue  terms  will 
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A,  r  <  a 


B,  r  >  a 


i  Fig  2.1 

Admissable  Deformed  Paths  of  Integration 
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occur,  except  one  corresponding  to  z  =  a  [since  this  is  a  simple  pole 
of  the  integrand].  Now  referring  to  Figures  2.1(A)  and  2.1(b)  we  make  the 
following  assumptions. 


(i)  T  represents  the  curve  A  a  a  a  A.  (as  before) 

(ii)  L  represents  the  curve  A  r  A  [in  particular,  we  choose 

the  straight  line  from  A  to  A  . ] 


(iii)  C  represents  the  curve  a  a  [the  indentation  around 
z  =  a  ] . 

X2 

,  as  before,  represents  the  integral  from  x^  to 

1 

[  ] 


x^  on  the  appropriate  curve  signified  in  [  ]  . 

Thus,  for  r  <  a 


Thus, 


ra 

pA  -i 

pA 

pa 

lim 

1 

+  / 

=  1-  + 

lim 

/ 

• 

r 

Ja 

■J  a 

-'A 

C 

J  a 

a  -»  a 

[r] 

[r] 

[L] 

a-*  a 

[C] 
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Now, 

fa  s(z)  (1-2rz+z2)6 

a  -*  a 

=  jti  g(a)  (l-2ar-KX2)6  , 

and  therefore, 


(2.5.3) 


K' 
2  jti 


— g(z) (l-2rz+z2)S  dz 
z-a 


g(a) (l-2ar+a2)S  + 


K' 

2jti 


1 

(z-a) 


g(z) (l-2rz+z2)8  dz 


For  r  >  a 


-2jti  (g(a)(l-2ar+a2)S)  , 


so  that 


2jti  (g(a)  (l-2ar+a2)S) 


Thus 


(2.5.4) 


K' 

2jti 


— g(z) (l-2rz+z2)S  dz 
z-a 


g(a)  (l-2ar+a  ) 


s 

+ 


K* 

2jti 


fK  T^y  8(z)(1"2rz+z2)s  dz 

[l] 
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Now,  by  Appendix  II 


2sS 


K'  g(a)  (l-2ar-Ki  )  =  1 


so  that  (2.5.3)  for  r  <  a  ,  becomes 


(2-5-5)  ifr  f[r]  Tifey  s(z)(i-2^«2)8  d* 


=  I  +  sfr  i^y  g(z)(1‘2rz+z2)S  dz 

[L] 


and,  (2.5.4)  for  r  >  a  ,  becomes 


(2-5-6)  sfr  ^r]  S(z)(l-2rz«2)s  dZ 


-  -  I  +  2^1  f-  s(z)(x-2rz«2)s  dz  - 


A 

[L] 


and  the  inversion  formula  (2.2.6)  (using  (2.4.6))  for  r  <  a  becomes 


(2.5.7)  H(r)  =  £  f  g(z)(l-2rz+z2)S  dz  , 


A 

[L] 


and,  for  r  >  a  becomes, 


(2.5.8)  H(r)  =  1  -  ~~  g(z)(l-2rz+z2)S  dz  , 


A 

[L] 


where  L  is  the  straight  line  from  A  to  A  . 
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§  6.  Change  In  the  Variable  of  Integration, 

The  integral 

(2.6.1.)  gfr  J_  JciSz)  g(z)(l-2rz+z"')S  dz 

A 

[L] 

where,  L  is  the  line  from  z  =  A  =  r-W  1-r2  to  z  =  A  =  r+is/l-r2  passing 
through  z  =  r  ,  has  a  critical  point  at  z  =  r  .  We  change  the  variable 
of  integration  so  that  the  critical  point  is  at  the  origin  in  terms  of 
the  new  variable  of  integration.  We  use  the  linear  transformation 


(2.6.2) 

z  =  r 

+  hJl 

-  r2  t 

so  that, 

(2.6.3) 

dz  = 

i si  1-r2 

dt  , 

(2.6.4) 

(a-z) 

=  a  - 

r  -  Wl-r2  t  , 

(2.6.5) 

(l-2rz+z2) 

=  (1  -  r2) (l  -  t2)  , 

(2.6.6) 

g(z) 

=  g(r  + 

1 

ro 

rt 

II 

rt 

V* 

and  L*  , 

the  t -image  of 

L  ,  is 

the  straight  line  from  t 

passing  through  t  =0  [i.e.,  a  segment  of  Im  t  =  0  ] . 

Thus,  the  t-integrals  corresponding  to  (2.5.7)  and  (2.5*8)  are 


for  r  <  a 
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(2.6.7)  H(r)  = 


K'  (1-r2)9*^ 

2k 


>+1 


-'-1  (a-r-Wl-r'' 

[L*  ] 


X(t)  (1-t2)8  dt 


and,  for  r  >  a 


(2.6.8)  H(r)  = 


K»(l-r2)S4^ 

2n 


,+l 


- 1  (  r  -a+i "V  1-r2  t ) 

[L*] 


x(t)  (l-t2)8 


dt 
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CHAPTER  III 

THE  ASYMPTOTIC  APPROXIMATION  FOR  DISTRIBUTION  FUNCTIONS  OF  TWO 
SAMPLE  SERIAL  CORRELATION  COEFFICIENTS 

In  this  chapter,  we  present  the  asymptotic  representations  of 
certain  integrals,  which  represent  distribution  functions  of  sample  serial 
correlation  coefficients.  The  method  used  in  obtaining  the  necessary 
asymptotic  expansions  is  the  same  as  that  which  is  presented  in  a  joint 
paper  by  A.  ERDELYI  and  M.  WYMAN  (1963)  entitled  "The  Method  of  Laplace." 
Notwithstanding  the  essential  similarity  in  the  developments,  there  is 
one  significant  modification  made  in  our  work.  Whereas,  in  Erdelyi  and 
Wyman  the  basic  functions  involved  in  the  Asymptotic  expansions  is  the 
Euler  integral  for  the  Gamma  function,  we  use,  as  basic  function  for  the 
expansion,  the  Euler  integral  for  the  Confluent  Hyper geometric  function. 


To  preserve  continuity  in  the  ensuing  presentation,  the  deriva¬ 


tion  of  the  complete  expansion  and  other  details  are  relegated  to  Appendix 
III.  We  shall  draw  upon  the  major  results  obtained  therein,  always  with 
appropriate  reference, 

§  1.  Basic  Assumptions  and  an  Asymptotic  Representation. 


The  integral  for  which  we  seek  an  asymptotic  representation  is 


[L*3 


where, 


'■'V.r z&id  &umh 


- 


- 

" 

:':oi  a;,  a  ;  'a  '  ■=>.:?:■  r  «rr  ,*•  •  £-* 
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(i)  L*  is  the  line  Im  t  =  0^ 

(ii)  s  is  a  real  variable  assuming  large  values., 

(iii)  a,r  are  real  variables,  a  fixed  and  such  that  |a|  <  1, 
and  r  varying  continuously  over  the  set  ((*:  |£|  <  1  and 

i  /  aK 

(iv)  x(0  is  an  analytic  function  of  t  regular  in  some 

neighbourhood  of  t  =  0  and  |  x( )  I  <  00  for  all  t  in 
[-l,+l]  i,e.  on  the  path  of  integration. 

By  Appendix  III, 

(3.i.2)  r+  — (i-t£)s  dt 

^  -1  a-r-i  s / l-r2t 

[L*] 


M  m 


a-r 

l-r‘ 


4-,1>  r<-^) 


nH-k-i  k 


s"  U  f  m+k+-| }  m+k+| ;  ^ 


m=o  k=o 


M  m-1 


X  7  ^  r(m+k+i)  {(T-ra2}  2  sU  U  (m-k+ijm+k+i;2^^2 


m=l  k=o 


+  0 


s  ( a~  r ) 2 

ll'rS 


) 


where,  U(  ;  ;  )  is  the  Confluent  Hyper geometric  function  in  the  notation 
of  SLATER  (i960),  and,  the  ckm>1^  s  are  given  by 

00  M 

(3.1.3)  y  y  ckTn’1^  sk  t2^1^  =  tx(t)+x(-t)]  exp[st2+slog(  1-t2)  ] 

m=o  k=o 


'Ar- 


{*)?■  ^-0 


;  v ,  i .  c 
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and , 


oo  m- 1 

(3.1.4)  ^  y  c<™  l’2^  sk  t2^nl+k  -  — p — -  [x(t)-x(-t)]exp[sts+slog(l-t'°)] 

m=l  to  2 


Now  by  equations  (2.6.7)  ancl  (2.6.8),  we  have,  for  r  <  a  , 


(3.1.5) 


H(r) 


M  m 


SL4rfL  t  t  ckm,1)  r<m+k+?)  -fc^T  sk 


m=o  k=o 


1  M  m- 1 
K'(l-r^S+2  y  V 

2jt  L  L 

m=l  k=o 


t1,2)  r(«k«) 


1  m+k- 


2  _k  „  f  ,  1  ,1 . s(a-r )2 

s  U  (  m+k+g » m+k+2 ; g 1 


+  0  ( fe^}M+2r(M+  i> u  ("+|  j  M+l  5  2fe§i~)). 


and  for  r  >  a  , 


(3.1.6)  H(r) 


=  1 


K 


M  m 

1  ( l-r2)S  Y  Y  _(m,l)  r ./  1  i\  f(r-a)2(m+  k  „  f  ,  1  s(r-a) 

2n  /  /  Ck  r(m+k+2)  |V^_r2  j"  s  U  m+k+  2 ; m+k+  2 ;  ^_rP 


m=o  k=o 


K1 ( l-r2)S+2 

2n 


+  0 


M  m- 1  1 

i  i  ym-i>2)  r(nH-k+i)  f&r " 


m+k+-| ;  tn+k+-| ; 


•a. 


1-r2 


m=l  k=o 


C  i> 
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Now, 


(5.1.7)  .  f 


and 


<7-  i 


{W+4 


exp[-sv]dv  , 


(3.1.8)  o(ii 


s(a-r)2 

exP  ~i-'r2 


/' 


v  2  exp[-v]dv 


ilaul 

l-r2 


Using  (3.I.7)  for  m+k  >  0  and  (3. 1.8)  for  m+k  =  0,  we  get,  in  terms  of 
the  above  integrals,  for  r  <  a  , 


(3.1.9) 


H(r) 


K'(l-r2)S  (0,1)  r/n  fs(a-r)21  P™  -k  r  n. 

=  — gj  o  ^2)  expl  V  l-r2~_,  /  exp[- v]dv 


s(a-r),2 

l-r2 


_i  Mm 

K1  ( l-r2)S~2  (q-r)  V  V  >,l)  „k 

2it 


ivi  III 

I  l 

m=l  k=o 


L  u  o 


00  m+k-i 

v _ _ 

5T 


exp[-sv]dv 


M  m-1 


K'(l-rg)S+S  V  y  (m-1,2)  k  T  P” 

2*  ^  7.  k  [  J  J(a- rjg  T 

m=l  k=o  °  l  l-r2  +  V| 


exp[-sv]av 


+  0 


M+i 

1  j  <- 


exp[ -sv]dv 


and,  for 


r  >  a  , 
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(3.1.10) 


H(r) 


K*  ( l-r2)S  (0,1)  n /  i  \ 

1  "  —  2^ -  Co  r^)  exP 


l(i-g): 

1-r2 


LjLr.qr.lf 

1-r2 


v  2  exp[-v]<lv 


K'U-r2)8'^  (r-g) 

2* 


M  m 

I  I 

m= 1  k=o 


(m,l)  k 


noo  m+k— 5- 

IjSFTJ 


exp[ -s v]dv 


K1 ( l-r2)S+2 

2it 


M  m-1 

I  E 

m= 1  k=o 


(m-1, 2)  k 


co  m+k-4 


°  4 


exp[-sv]dv 


+  0 


M+i 


°  {-W +  4 


exp[-sv]dv 


. 


Thus  as  r  ->  a  (3. 1.9)  becomes 


(3.1.11)  H(a-O)  =  c^0’1) 


K^(l-g2)S+5  Y  V  c*(m- 1,2)  r(nvfk-2  ) 

2jt  /  ,  /  k  m-A 


M  m-1 


m=l  k=o 


+  0,73. 


and  (3.I.IO)  becomes 
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(5.1.12)  H(a+0)  .  1  -  c^0’1) 


i  M  m- 1 

K*(l-g2)S+2  y  V  *(m-l,2)  r(nn-k-j) 

2k  /  ,  /  ,  k 


m=l  k=o 


+  0 


M+|r 


where, 

(3.1.13) 


K*  =  K'  at  r  =  a  . 


(5.1.14) 


:*(  '  '  =  cf  ’  '  at  r 

k  k 


=  a 


and  the  interchange  and  limits  and  integration  is  justified  since  the 
integrals  are  uniformly  convergent. 

Now  by  (3.1.3) 

(3.1.15)  .  x*(0)  . 

where,  by  equation  (2.6.6) 


X*(0)  =  g(a) 


and 


(3.1.16)  x(fc)  =  X*(t)  at  r  =  a  . 


By  Appendix  II,  (equation  15) 


K*(l-q2)a  X*(p)  _  K*(l-qa)Sg(q^  1 

2  2  =  2 


(3.1.17) 


(  ,  )  ii 

'  Si.  1  i  :  <  •  .  : . 

■  'i  s'  ,,  f 

f,  1 


)l  I. 


Thu.; , 


(3.1.18)  H(a-O)  =  H ( a+O )  =  H(a) 


K*( l-a5)s+^ 

2it 


I 


ni=l 


m-  1 

Y  c*(m-l,2) 

k=o 


r(rn+k-^) 

s  2 


+  0 


We  take  (3. 1.18)  to  be  the  asymptotic  value  of  H(r)  at  r  =  a.  It  may 
be  noted  that  the  expansion  (3. 1.18)  may  also  be  derived  independently, 
directly  from  the  integral,  by  the  method  outlined  in  ERDELYI  and  WYMAN 

(1963). 

Thus  in  equations  (3*1.9)  and  (3.1.10),  for  M  =  0,  we  have 
the  following  asymptotic  expressions,  for  r  <  a 


(3.1.19)  H(r)  =  S'.iA^JS  x(0)  r(4)  J 


1-r2 


v  2  exp[-v]dv 


+  0 


1 

v2 


fei2  +  4 


exp[-sv]dv 


:or  r  >  a 


(3.1.20)  H(r)  =  1  -  x(0)  r(4)  expr^|^|  j"  v"2  exp[-v]dv 


1-r2 


+  0 


(L 


f(  r-af  "I 

0  1i^-  +  4 


exp[ -sv ]dv 


which,  for  r  =  a,  give 
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(5.1.21) 


which  would  be  an  approximation  to  the  first  term  in  (3. 1.18),  Thus,  we 
take  as  the  asymptotic  representations  for  r  <  a  , 


(3.1.22)  ll(r)  ~  x(o)  r(i)  exp 


and,  for  r  >  a  , 


s(q-r )2 

'  l-rs 


p°°  __i 

I  v  2  exp[-v]dv 

s(q-r )g 

1-r2 


(3.1.23)  H(r)  ^  1  -  K'  x(°)  r(J)  expp^"^  P  v  2  exp[-v]dv  . 


kis£ 

1-r2 


By  equation  (2.6.6), 


(3.1.24) 


x(o)  =  g(r)  • 


Thus  corresponding  to  equations  (2.5*7)  and  (2.5*8),  we  have  for  r  <  a 


(3.1.25)  H(r)  ^  l^hL^Ljdil  r(^)  exp[lki|lf  j'  v"2  exp[-v]dv  , 


(a-r  V 

1-r2 


and  for  r  >  a. 


(3.1.26)  H(r)  1  -  -  r(i)  expF-^-^-  J  v"2  exp[-v]dv 


Liz <£. 

1-r2 


, ...  «  o 


F 
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§  2,  Approximations  to  Distribution  Functions  of  two  Sample  Serial 

Correlation  Coefficients. 

By  equations  (2,^,6),  (2.5.7)  an<*  (2.5*8),  the  z-integral  has  form 

(3.2.1)  |jl  j r  ~  g(z)  ( l-2rz+zs)s  dz 

A 

for  both  the  circular  and  non-circular  cases. 


(i)  The  Circular  Case. 


By  equation  (2.3.20),  the  correspondence  with  (3.2. l)  is 


(3.2.2)  (1)  s  =  |  -  1 

(il)  K'  =  - - - 

R 

( l~2ar+a2)2 

(111)  g(z)  =  J— 


Thus  by  equations  (3. 1.25)  and  (3. 1.26),  for  r  <  a 


(3*2.3)  H(r)  ^ 


2  s fit: 


(1-v2) 


exp[(|  -1) 


iocll 


r2  -1 


( l-2ar+a2) 


l-1 


( 1-ar ) 


and,  for  r  >  a 


-  Ii4 


(3.2. U)  H(r)  1  -  -L~ 

2  si  jt 


(l-r.2J 


( l-2ar+a2)' 


-1 


exp 


U  _D  Lira)-. 
L >2  ;  1-r2  J 


( 1-ar) 


and,  in  particular,  for  r  =  a  , 


{(1  -D^} 


exp[-v]dv 


(3.2.5) 


H(a)  a.  |  . 


( ii^  The  Non-circular  Case. 


By  equation  (2.3.21),  the  correspondence  with  (3.2. l)  is 


(3.2.6) 


(i)  8  =  |  -  1 


(n)  k'  =  — 


( l-2ar+a2) 


S  -2 
2 


2 

2  2 

(iii)  g(z)  -  ( i-az)5( 1-ar+az-rz) 


Thus,  by  equations  (3.1.25)  and  (3. 1.26),  for  r  <  a  , 


(3.2.7)  H(r)  -v  — (.1^2)^  Ll~rsl 


n  1 
2  "  2 


2 


(1-ar)2  (l-2ar+a2) 


—  -2 

2 


exp 


(f  -1) 


r 


{q  -D^ 


v  2  exp[-v]dv 


x 


ft,  ,  ftrfti 
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and  for  r  >  a 


(3-2.8) 


H(r)  -v  1 


n  1 

_J_  (l-a2)*  (l-rg)2'2 

2  v/rt  §  -2 

( 1-ar  )2(  l-2ar+a2) 


x 


exp[-v]dv 


9 


and,  in  particular,  for  r  =  a  , 

(3.2.9)  H(a)  ~  | 


The  magnitude  of  error  involved  in  making  the  above  approximations, 
for  both  the  circular  and  non-circular  cases,  is 


(3.2.10) 


This  order  relation  is  uniform  in  r  and  n,  A  more  exact  relation 
uniform  in  n  and  not  in  r  may  be  derived  from  the  equations  (3,1,19) 
and  (3.1.20). 


Now,  for  r  <  a,  consider 


(3.2.11) 


exp 


(f  -1) 


(a-rV 


/' 

{<f  -D1! 


v  2  exp[-v]dv 


(a-r); 


Put 


•  ■  ' 


c  3  :  i ",  • 


;l.  "  ‘  .  .  ■ 
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Then , 

dv  =  t  dt  , 

vi/2/t 

and, 

v  =  (“  -l)  ^v~r-^  -  corresponds  to  t  =  N^n~-  (a— Ll 
2  x-r 


where  $(  )  is  the  standardized  normal  distribution  function  and  cp(  ) 
the  standard  normal  probability  density  function. 

The  corresponding  result  for  r  >  a  is 


r-  (r  -nr  ^ 


L  \l  1-r2  -1 


(3.2.13) 


n.SiiT 
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Thus,  for  the  circular  case,  we  have  by  equations  (3.2.3)  and  (3*2.4), 
for  r  <  a  , 


(3.2.14)  H(r) 

n/2* 


ilirfl 


5  -1 


( 1-ar )( l-2ar+a2)2  ^ 


-  -J n-2  la-r| 

\Tl-rg 

-  v/n-2  la-r  1 


and,  for  r  >  a 


(3.2.15)  H(r)  ~  1  - 

v  2it 


XL-rf) 


••  J  n-2  la-r| 


( l-ar )( 1  -2ar-t-a2) 


S-i 
2  cp 


-  \/n-2  |a-rl 

.  sJT^P 


For  the  non-circular  case,  we  have  by  equations  (3.2.7)  and 
(3.2.8),  for  r  <  a 


1  -  -  A 

(3.2.16)  H(r)^—  — OftF  (■1_r2) . . 

§  -2 

( 1-ar  )2(  l-2ar+a2)2 


\/n-2  !a-r| 


vT 


cp 


sfn-2  la-r  | 


and,  for  r  >  a 


(3.2.17)  H(r)  ^  1 


n  _  _1  ^ 

1_  (l-a2)^(l-r2)2  2 

n-2 

( 1-ar )2( l-2ar+a2)2  ^ 


-  \fn-2  la-r | 

'  sfwP 


n/s rt 


-  \/n-2  la-r | 


r 
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As  discussed  earlier,  in  equation  (3.2. 10) ,  the  order  of  the 
magnitude  of  the  error  involved  in  making  the  above  approximations  is 

(3.2.18)  R  =  0  ^  n‘2 

uniformly  in  n  and  r. 

These  results,  namely,  (3.2.14),  (3.2.15),  (3.2.16)  and  (3.2. I7), 
are  represented  graphically  in  Appendix  four.  The  order  of  the  magnitude 
of  error  has  not  been  taken  into  account  in  the  tracings  of  the  curves. 
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APPENDIX  I 

COMPLEX  INTEGRALS  IN  THE  SENSE  OF  A  CAUCHY  PRINCIPLE  VALUE 


Assumptions : 

(1)  f(z)  is  a  function  defined  in  some  region  of  the  z-plane. 

(2)  Consider  a  fixed  point  z  =  a  ,  and  for  some  R  >  0  (fixed), 

consider  a  neighbourhood  |  z-cx- 1  <  R.  f(z)  is  taken  to  be  regular  for 
all  points  of  |z-a|  <  R  except  at  z  =  a. 

(3)  The  point  z  =  a  is  a  simple  pole  of  f(z)  . 

(4)  C  is  any  continuous  curve  in  the  z-plane  passing  through  s  =  a  and 

intersecting  with  |z-a|  =  R  at  z  =  A  and  z  =  B. 

We  wish  to  determine  sufficient  conditions  under  which  the  integral 


lim 
a  -*•  a 

b  ->  a 


f(z)  dz 


(1) 


exists  in  the  sense  of  a  Cauchy  Principle  Value,  where  the  limits 
a  — >  a,  b  a  are  taken  to  be  limits  through  the  set  C. 

The  configuration  is  as  in  the  diagram. 
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By  assumption 


f(z)  =  ~T  +  g(z) 
z-a 


for  all  z  in  |z-a|  <  R  where  g(z)  is  regular  for  all  z  in 
| z-a|  <  R,  and  K  is  constant  with  respect  to  z.  Thus, 


lim 
a  -*■  a 
b  — ►  a 


r 


j 


f(z)  dz 


lim 
a  -+■  a 
b  -v  a 


B\  r 


+  g(z) 
z-a  ' 


dz 


.  u»  ( r\  rBv Ldz+  lim  (r\  r 

!\  k.  «  kJ  a  '  z— a  a  v  ,v  ^  A.  ^  b 


a  — *■  a 

b  ->•  a 


a  -*•  a 
b  ->  a 


g(z)dz 


Now,  by  regularity  of  g(z)  for  all  z  in  |z-a|  <  R  , 


and,  further, 


lim 
a  -»•  a 
b  — ►  a 


r 

^  A 


^  ^  g(z)ds 

J  b  7 


rB 

=  /  g( z)dz  , 

J  k 


lim 
a  — ►  a 
b  ->  a 


c  r  *  r ) 


K 

z-a 


dz 


=  K  lim 
a  — >•  a 


log 


fa -a  e  B-al 
[A-a  b-a| 


K  lim  log  ( -r—  }  +  K  lim  (r“" 

a  —v a  \A~a/  a  — ►a  \b-<X/ 

b  a  b  — ►  a 


K  log  (~)  +  K  lim 

V A-a/  a  -►a  \b-a, 

b  -►  a 
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Thus 


b  -»■  a 


where  the  integral  on  the  left  will  exist  if  and  only  if 


a  -v  a 
b  — ►  a 


is  finite. 


This  is  a  sufficient  condition  that  (l)  exist  as  a  Cauchy  Principle  Value. 
This  above  condition  may  be  transformed  to  something  much  simpler  in  form. 


Consider  the  normal  to  the  curve  C  at  z  =  a,  and  consider 
the  transformation 


z  -  a 


t  e 


id 


where  z,t  and  a  are  complex  variables,  9  is  a  fixed  angle  determined 
so  that  the  normal  to  the  curve  C  at  z  =  a  in  the  z-plane  is  trans¬ 
formed  to  be  the  line  Im  t  =  0  in  the  t-plane.  Thus  in  the  t-plane, 
the  limit  becomes 

t 

lim  K  log  ~ 

t  —0  b 

a 

t,  -»0 
b 

where  t^,  t  are  the  t-images  of  a  and  b  respectively,  and  the 
limits  are  through  the  set  C  (the  t-image  of  C  in  the  z-plane). 


v  /;■-  :  .  o 
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Now  we  consider  a  special  case  of  interest.  Assume 
the  t-image  of  the  curve  C  in  the  z-plane  is  symmetric  about 


and  that  t  =  t, 
a  b 


(the  complex  conjugate  of  t^ ) , 


Thus  if 


fcb  =  iv  6 


t  =  t,  e 
a  1  b  1 


icp 


“icp 


0  <  cp  <  2it 


and 


lim 

t  ->  0 

a 

tb  -0 


becomes 


lim 

cp  -►  <5 


Thus  we  have 


lim  K  log  ( t  /t,  j  =  lim  K  log  fe 
t  — ►  0  a  '  cp  -»•  $ 


t,  — ►  0 
b 


=  lim  K  log  (  e 

cp  — >  $ 


~2icp^j 


=  lim  -2iKcp 
cp  -»  $ 

=  -2iK<3> 


which  is  constant  and  finite. 


[cf .  SANSONE  and  GERRETSEN  ( i960)  pp .  127-128], 


that 

Im  t  =  0, 
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APPENDIX  II 

CERTAIN  EQUALITIES 


Equations  (2.3.18)  and  (2.3.19)  are  of  the  form 


(1) 


cp(C(z))  C  (z) 

£(z) 


h(z)  (l-2rz+z2)S 

z-a 


where  s  is  defined  as  in  equation  (2.4.3). 


In  equation  (2.3. 18) 


K 


( l-2ar+a2)2 


-1 


and 


h(z) 


_ 1-z2 

( l-az)(  l-zn) 


y 


while  in  equation  (2.3,19)» 


,  --  d-a2)2 

n 

( l-2ar+a2)2 


and 


h(z) 


2 

— 1L-Z fll 

( l-az){  ( l-az)2(  l»ar+az-rz)2+z2n  2(a“z)2(a-r+z-arz)2}2 


and  in  both  equations  (2.3. 18)  and  (2.3. 19) 

(2)  L'±z).  =  _( I ' -zfl.il 

'  £  (z)  ( z-a) ( 1-az) ( l-2rz+z2) 
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Thus , 

(5)  Kh(!)(iw)'  =  *(»*» 

and, 

(4)  K  h(a)  (l-2ar+a2)S  =  lim  K  h(z)  (l-2rz+z2)S 

z  — ►  a 


-  ^(^(a))  [  i-a2  1- 


-2ar±af 

2ar+a2 


}  =  <P«(a))  . 


Now, 

(5) 


£(z) 


iz( l-2ar+a2) 
l-2rz+z2  * 


so  that 


(6) 


£(a) 


ia(  l-2ax+a2) 

l-2ar+a2 


=  ia  -  ia  =  0 


Thus , 

(7)  K  h(a)  (l-2ar+a2)S  =  q>( 0)  , 

and  since  Cp(  )  is  a  characteristic  function 


(8)  cp(0)  =  1. 

Thus, 


K  h(a)  (l“2ar+a2)S 


(9) 


=  1  . 
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Comparing  (2.5.20)  with  (2. It. 5)  we  see  that 


(10)  K'  g(z)  ( l-2rz+z2)S  = 


( l-2qr+q2) 


I-1 


1-z2 

1-az 


( l-2rz+z2)' 


1-  z 
l-ar 


1-a 


l-z‘ 


( l-2ar+a2)‘ 


5-1  (l-az)(l-zn) 


( l-2rz+z2)‘ 


1-z 

l-ar 


K  h(z)  (l-2rz+z2)S 


Thus, 


(11) 


K*  g(a)  (l-2ar+a2)S  =  1  ~~~  K  h(a)  (l-2ar+q2)S 


1-a 


n 


=  1  . 


Again,  comparing  (2.3.21)  with  (2.4.3),  we  have 

l 

(12)  K'  g(z)  (l-2rz+z2)S  =  ■  . 


2 


( l-2ar+a2)' 


^  -2  ( l-az)2( 1-ar+az-rz) 


( l-2rz+z2) 


■A 


f  ( l-az)2(  l-ar+qz-rz)2+z2n  2(q-z)2(q-r-)-z-arz)2}; 

( 1-qz )  (l-ar+az-rz) 

I  2 

( 1-a2)2  _  (1-z2)2  (l-2rz+z2)S 

H  «°2  2n —2  <-  JL 

(l-2ar+q2)2  ( l-az){  ( l-qz)2(  l”qr+qz“rz)2+z  (a-z)2(a-r+z-arz)2]2 
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f  ( l-gz)3(  1-gr-NXZ-r z )2+z^n  ^q-z^g-r+z-grz)2]^ 

(l-gz)  (1-ar+az-rz) 


K  h(z)  (l-2rz+z2)S 


Thus 


(13)  K'  g(a)  (l-2ar+a2)S  =  l'-'far^}  K  h(a)  (1‘2ar+a'2) 


=  1  , 


which  justifies  the  statement  of  equation  (2.5.4(a)).  Thus,  from 


(14)  K1  g(a)  (l-2ar+a2)S  =  1 


we  conclude  that 


(15)  lim  K1  g(a)  (l-2ax+a2)S  =  K*  g*(a)  (l-a2)S 


where 


K'  =  K*  for  r  =  a 


and 


g(z)  = 


g*(z)  for  r  =  a 
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APPENDIX  III 


AN  ASYMPTOTIC  EXPANSION. 


In  this  appendix  we  give  the  details  of  the  derivation  of  the 
asymptotic  expansion  for  the  integral 


(1) 


J 


.+1 


X(t) 


(l-t2)8  dt  , 


-1  a-r-Wl-r2  t 

[L*] 


with  the  assumptions  (3.1.1)  (i)  -  (iv). 


Thus  consider 


(2) 


I(r)  = 


C 

[L*] 


Alii 


a-r-Wl-r2  t 


(l-t2)S  dt 


n/T 


1  r+] 

:p~  J-i 


[L*] 


x(t) 


it 


(i-t2)s  at 


p  +1 

J-  1 

[L*] 


x(t) 


a~r 
si  i-r2 


+  it 


{ 4 

l  l-r2  J 


(l-t2)S  dt 


r 


^+i 


si  l-r2  L  -'o 

[L*] 


x(0  (  7== 

U fu? 


.+  it 


(a-r  { 

l-r2 


+  t 


(l-t2)S  dt 


+1  x(“°^  s/T^"  1C 


J  o 

[L*] 


;is«f 

l-r 


+  t‘ 


(l-t2)S  dt 


■ 
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a- 


L  >/T-] 


o 

[L*] 


U(0  H  x(-03 

-r  )2 


f 

L  l-i 


+  t‘ 


(l-t2)S  dt 


+ 


^(0  -  x(-Q]  (1.t2)s  dt 


f 


each  integral  being  absolutely  convergent. 

Referring  to  equation  (2.6.6),  we  have, 


(5)  g(z)  =  g(rWl-r2  t)  =  x(t)  , 

where  g(z)  is  an  algebraic  function  of  z  ,  which  is  real  valued  for  all 
real  values  of  z  ,  and  is  regular  in  some  neighbourhood  of  z  =  r  .  Thus 
for  some  R  >  0  and  |z-r  |  <  \l  1-r2  R 

00  00 
(4)  g(z)  =  \  a 

m=o  m=o 


where  the  a  are  real.  Thus, 
m  7 


00 

(5)  X.(t)  +  x(-t)  =  2  (-l)m  a2m  (l-r2)m  t2m 

m=o 


and 


J 


(6) 


00 


Put 

(7) 

and 

(8) 

which 

(9) 

and 

(10) 

where 

(u) 


t(x(t)  -  x(-t)l  =  21  y  (-1)'"  a2m+1  (l-r2)m+*  t2m+2 

m=o 


=  ^  (-l)™  a2m+1  (1-r2)"’  t2m 

m=o 


x(0  +  x(-t)  =  2  f^t2)  , 


t[x(t)  -  x(-t ) ]  =  2vJ 1-r2  t2  f2  (t2) 


means  that  for  |t|  <  R 


f^t2)  =  £  (-l)m  a£m  (l-r2)m  (c2)m 

m=o 


f2(t2)  =  £  (-1)"  a2m+1  (l-r2f  (t2)m 

m=o 


a 

m 


is  always  real. 


Thus  (2)  becomes 
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(12)  I (r )  =  2fa‘rl 


1  -  r‘ 


J  o 
[L*] 


fx(t:2) 


{^*•4 


(1-t2)1  cl  L 


-  2 


J  o 

[L*] 


C£  £g(t2) 


14 


(i-t2)8  dt 


a-r 


1-r2  Jo 
[L*] 


Sax. £  +  t* 

l-r2 


(l-t8)s  d(t2) 


'  o 
[L*] 


t  f2(t2) 


+t. 


l-r- 


(l-t8)S  d(t2) 


a-r 

l-r2 


J  o 
[D] 


fa-tY 

1  -r2 


+  v 


(l-v)S  dv 


J  o 
CD] 


v2  f2(v) 


(a-r ); 

l-r2 


+  v 


(l-v)S  dv  , 


where, 

(13) 


v  =  t£ 


and  D  is  the  v-image  of  L* 


-  64 


Now  consider 

(14)  (l-v)S  =  exp[s  log(l-v)] 

=  exp[-sv  [  -  ] 

and  let 

(15)  .loslkii  .  f(v)  . 

Now 

°°  k 

(16)  f(v)  =  ,  (|v|  <  l)  . 

lc=o 


Thus  ^(v)  is  regular  at  v  =  0  ,  and  thus  is  regular  for  all  v  in 

[0,1)  . 

Thus  (12)  becomes, 


(IT)  l(r )  = 

1-r2 


where , 


(i)  f£(v)  ,  i 


=  1,2,  and  f(v)  are 


exp[  ~sv  V/(v)  1  dv 


exp[-sv  ip{y)  ]  dv  , 


regular  in  a  neighbourhood 


of  v  =  0 
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(ii)  Re(sv  V/(v))  >  0  for  all  v  in  (0 , 1  ]  .  [Indeed 

sv  ip(v)  is  real  valued  and  monotonically  increasing 
in  [0,1]  ,  which  is  a  path  of  "steepest  descent"], 

(iii)  for  some  fixed  /  0  ,  and  in  [0,1)  ,  there 

is  a  fixed  5  >  0  ,  such  that 

IV'(v)  |  >  B 

for  all  v  in  [vq,1]  ^ 

(iv)  the  integrals  are  absolutely  convergent. 

Thus  we  concern  ourselves  with  deriving  the  asymptotic  expansion  of  an 
integral  of  the  form 

(18)  I*  =  f  — - ~iyJ —  exp[-sv  ^/(v)]  dv 

(a-r)2  \ 

—~r  +  vf 

l-r2  J 

under  the  conditions  (17), (i)  -  (iv). 

From  ERDELYI  and  WYMAN  (1963)  we  have  the  following! 

Assume  that  the  integrand  has  the  form 

(19)  H  =  v^  ^  f*(v)  exp[ -sv^(v)  ] 
with  the  following  assumptions; 

(20)  (i)  The  v-plane  is  cut  along  arg  v  =  7  from  v  =  0  to 

iy 

v  =  00  e  ,  7  <  arg  v  <  7  +  2tc 
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(ii) 

p  is  a  fixed  real  number  and  p  >  0  . 

(iii) 

A  is  always  a  fixed,  finite,  complex  number.  When 

Re  A  >  0  ,  p  >  0  suffices,  when  A  unrestricted  it 

is  necessary  that  p  >  1  . 

(iv) 

f*(v)  and  ^/(v)  are  both  analytic  on  the  path  of 

integration,  say  D  ,  and  for  some  fixed  R  >  0  ,  are 

both  regular  in  |v|  <  R  . 

(v) 

s  is  a  complex  number  -*  oo  through  the  closed  unbounded 

point-set,  S(A)  ,  of  the  s-plane,  where  S(A)  is 

defined  by 

(4f-l)  |  -  py  +  A  <  Arg  s  <  (4f+ 1)  |  -  py  -  A 

where  A  is  fixed,  0  <  A  <  jt  . 

(vi) 

For  all  paths  of  integration,  D  ,  used  /  Hdv  exists 

as  an  absolutely  convergent  integral  for  each  sufficiently 

large  s  in  S(A)  „ 

Under  these  conditions  the  following  theorem  holds . 


Theorem  l; 

Let 

be  a  continuous  curve  joining  the  fixed,  finite,  point 

v  =  ^  0  to  the  fixed  point  v  =  B  ,  where  B  might  or  might  not  be 

the  point  at  infinity.  Further  suppose  that  the  points  of  D-^  are  contained 


in  the  point-set  |v|  >  |v^|  .  If 


(21)  H  =  vpA  ^  f*(v)  exp[ -sv^(v)  ] 

satisfies 


(i)  |^(v) |  >  5^  >  0  ,  for  some  fixed  5^  ,  uniformly  for  points 

of  . 

(ii)  for  t(s,v)  =  Arg(svP  ^/(v))  ,  cos  t  >  £>2  >  0  ,  for  some 

fixed  >  uniformly  in  (s,v)  ,  where  s  is  sufficiently 

large  in  S(A)  and  v  belongs  to  . 

(iii)  /  Hdv  exists  as  an  absolutely  convergent  integral 

jl\\ 

for  each,  fixed,  sufficiently  large  s  in  S(a)  , 


then 

(22) 


0 


uniformly  in  s  ,  as  s  «-»  00  in  S(A)  . 


The  order  relation  implied  in  (22) 


1  .e . 


[D] 


Hdv  =0 


f  1 


s  p 


for  n  =  0 ,  1,  ... 


is  uniform  in  s  only,  since  there  is  no  assumption 


that  the  uniformity  extends  to  any  parameters  that  may  occur  in  the  integral, 


[  cf  .  ERDELYI  and  WYMAN  (  19 63)  and  (1963a  ) 
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In  this  appendix,  the.  following  restrictions  are  maintained  with 
reference  to  20, (i)-  (vi)  . 


(23) 


(i)  the  v-plane  is  cut  along  arg  v  =  0 


(ii)  p  =  1 

(iii)  A  is  real  and  always  A  >  0 

(iv)  the  same  as  in  (20) 

(v)  s  is  a  real  variable  — »  oo 

(vi)  all  paths  of  integration  are  along  the  line  Int  v  =  0  , 
from  v  =  0  to  v  -  B  ,  inhere  B  might  or  might  not  be 
the  point  at  infinity. 


Digression  on  the  Confluent  Hyper geometric  function. 


We  adopt  the  notation  of  SLATER  (i960).  All  variables  and  parameters 
(except  perhaps  the  variable  of  integration)  are  assumed  real.  Thus  for 
a  >  0  ,  x  >  0 


1-b  p  00  a- 1 

(2k)  U(a  ;  b  ;  x)  -  f^y  J °  dt  • 


Set  x  =  sA  ,  where  s  >  0  is  real.  Thuss 


(25) 


u(a  ;  b  ;sA)  = 


1-b  . 1-b 
s  A 

r(a) 


a-1 


(sA+t } 


a-b+1 


exp[-t]  dt  . 


Set  t  =  sv  ,  s  >  0  as  above  . 


U(a  ;  b  ;  sA) 


exp[-Gv]  dv 


Thus , 
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(26) 


1  -b  n  oo  a- 1 

^  /  v 


f(a)  -'o  {A+v} 


a-b+1 


For  our  purposes, 

(i)  b  =  a  ,  and 

(ii)  A  is  a  real  valued  function  of  arbitrary  parameters, 
such  that  A  >  0  . 

Thus,  for  a>0,  A  >  0  ,  we  have 


(27)  U(a  ;  a  ;  sA)  = 


,1-a  poo  a-1 

^  /  v 


r(a)  -'o  {A  +  v} 


exp[-sv]  dv  , 


or 


(28) 


P°°  a"1  .1 

/  — -  exp[-sv]  dv  =  r(a)  Aa  U (a  *  a  *  sA) 

o  {A+v} 


Now  consider  the  behaviour  of  U(a  *  a  $  sA)  as  s  ~>  °o  .  It  can  be  show 


m 


that 


(29)  u(a  ;  a  ;  sA)  =  -1—  {1  +  0(~)) 


a. a  v  '•sA/ 

s  A 


and  thus. 


(30) 


a-1 

f  — 

-  0  [A+v } 


exp[-sv]  dv  = 


pi  a  +°<i»  • 

A  s 
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Application  of  Theorem  I  to  U(^A+q  ;  A+q  \  -s-^a~r )  \  q  =  0,1,... 

\  1-r2  / 


We  consider  the  integral 


(31) 


A+q-1 


J  o 
[D] 


exp[-sv]  dv  = 


=  T(A+q)  | 


III 


2) A+q-1 


v(Mq  ,  X+q  ;  iiSCl£ 

\  1-r2 


For  any  fixed  v  ^  0  on  D  lying  between  0  and  °o  ,  J  +  vi- 

°  l  1-r2  J 

is  uniformly  bounded  away  from  zero  for  all  v  e  [v  ,  oo)  i.e.,  uniform  in 

the  parameter  r  .  [We  recall  that  r  varies  continuously  over  the  set 

(I-  :  |£  |  <  1  and  i-  ^  a)  ,  where  a  is  fixed.]  Thus  -f  +  v|-  is 

l  1-r2  ' 

uniformly  bounded  for  all  v  in  [v  ,  oo)  .  Thus  the  integral 


(32) 


/' 

J  n 


A+q-l 


exp[-sv]  dv 


with  the  identification 


f*(v) 


{ 


(cL-r  y 

1-r2 


+  v 


-1 


and 


Tp(v)  =  1 


satisfies  all  the  conditions  of  Theorem  I  . 


3 


Thus, 


(33) 


o 

[D] 


v 


A+q  - 1 


(g-O1 

l-r2 


+  v 


exp[-sv]  dv  £  0 


and 


rc 

(34)  / 


CD] 


v 


A+q-1 


(Sizll 

l-r2 


+  V 


exp[-sv]  dv  £ 


[D] 


A+q-1 


(a-rV 

l-r2 


+  v 


exp[-sv]  dVj 


uniformly  in  s  and  r  ,  by  virtue  of  the  uniform  boundedness  of 


(g-r)- 

l-r2 


-1 


+  vr  for  all  v  in  [v  ,  . 


Application  of  Theorem  I  to  the  integral  I*. 


We  consider  the  integral 


(35) 


I*  = 


V*'1  thl 


J  o 
[D] 


(g-r); 

l-r2 


+  v 


exp[-s v^(v)]  dv  of  equation  (l8) 


which  satisfies  conditions  (17),  (i)  -  (iv)  . 


With  the  identification 


f*(v)  = 


(q-r  V 

l-r2 


+  V 
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I*  satisfies  the  conditions  of  Theorem  I 


Thus, 


(56) 


o 

[D] 


A-l  w  ^ 
v  f(v) 


(a-r); 

l-r2 


+  V 


exp[~sv  V/(v)]dv  %  0 


and 


(57) 


I*  = 


A-i  . 
V  f(v) 


J  o 

[D] 


(a~r)‘ 

l-r2 


+  v 


exp[ -sv  ^(v) ]  dv 


J  o 

[D] 


A-l 


v 


=El 


— (_l)_ exp[-sv  ^/(v)]  dv  , 


+  v 


uniformly  in  s  and  r  by  virtue  of  the  uniform  boundedness  of 

-1 


Us=iE  + 

L  l-r2 


v 


for  all  v  in  [v  ,  1] 

o’ 


The  Asymptotic  Expansion  for  the  Integral  I*. 


We  consider  the  integral 


(58) 


rV 
/  ° 
O 

:d] 


A-l  v 

V  f(v) 


f 


l»l  l  1. 


+  V 


exp [-sv  ^(v)]  dv 


where  we  choose  vq  fixed  such  that  f(v)  and  ^'(v)  are  regular  for 
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I v |  <  2 | vq |  <  1  .  Write 


(39)  exp[-sv  ip(v)  ]  =  exp[-sv]  exp[sv  (l-^(v))  ] 


Then, 


(1*0) 


A-l 


V  ~  f  ( V  ) 


• 0  f  (°-r ) 

[Dl  l  1-r2 


+  V 


exp [ -sv  V'(v)]  dv 


V°  f(v) 


[D] 


1  1-r 


+  v 


exp[sv  {l-^/(v)}  ]  exp[-sv]  dv 


iV 


o 

[D] 


*-l  -/  \ 

V  f(v) 

{<5?  *4 


G(sv,  v)  exp[-sv]  dv  , 


where, 

(41)  G(sv, v)  =  f(v)  exp[sv  (1  -  ^(v)}  ] 

Now  put  sv  =  w  and  treat  G(w,  v)  as  a  function  of  two  variables.  Thus 
for  | v |  <  2 | vq | 


f(v) 


L  b-v’ 

n=o 


and, 


(^3) 


r.f(v)  .  1+l£iikzl  .  -  £  £ 

n=l 


and , 


(44) 


which 


(**■5 ) 


where 


(46) 


where 


(^7) 


Since 
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exp[w  (l-V'(v)) 


y. 

n=n 


n! 


00  00 


on  rearrangement,  gives  the  pox^er  series 


^  Qn(w)  vn  , 

n=o 


Qn(w) 


is  a  polynomial  in  w  of  degree  <  n 


Thus, 


f(v)  exp[w[l-^(v))] 


00 


n=o 


00 


{I  b" v"}  { I 


Qn(w) 


n 


v 


n=o 


=  X  p»(w) 

■n=o 


LI 

pn^  -  X  b"-k  Qk(w) 

k=o 


n 


■  1 


(n)  k 
c,  '  W 
k 


k=o 


for  I v I  <  2|vo|  ,  all  component  series  converge 


,  ji-fetg] 

n=o 


P  (w)  vn 
nv  ' 


:  >  r '  ?:> 
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converges  to  f(v)  cxn[w ( l-^(v) ) ]  and  is  the  Taylor  series  expansion  of 
G(w,v)  . 


Thus  the  finite  Taylor  expansion  is 


(48) 


G(w,v) 


n=o 


|v|  <  |v 
i  i  _  i  o 


where, 


V. 


N 


fN(w»v) 


N+l 

v 


t 


and, 

00 

fN(w>v)  =  ^  (W )  vn'N_1  . 

n=N+l 

Now,  G(w,v)  is  uniformly  bounded  and  regular  in  v  on  D  for  |v|  <  2jvQ 
and  all  s  .  Thus  f  is  a  uniformly  bounded  and  regular  function  of  v 
for  | v|  <  2 | v^ |  on  D  . 

By  means  of  the  Cauchy  inequalities  the  magnitude  of  the  remainder 
V^T  may  also  be  expressed  by  means  of  an  order  relation. 

Since  ^/(v)  is  regular  for  |v|  <  |v  |  and  ty(0)  =  1 
|l-^(v) |  <  K|v  |  ,  for  some  fixed  K  >  0  ,  uniformly  in  v  for  |v|  <  2  j  ] 
providing  |v  |  is  taken  to  be  fixed  and  sufficiently  small.  Thus 
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m 


=  0  (vN+1  exp[l(|vo|  |w|])  ns  v  -*  0  ^ 


Thus,  returning  to  (40)  we  have 


(50) 


J  o 
[D] 


A-l 


(q-r)> 

l-r2 


+  v 


G(sv,v)  exp[-sv]  dv 


N 


■  l  I 


n=o 


o 

[D] 


P  (sv)  V 
n'  ' 


A+n-1 


fiscjOf  + 

l  l-r2 


exp[-sv]  dv 


P 

J  o 

[D] 


A-l 


v 


(q-r); 

l-r2 


+  v 


exp[-sv]  dv 


© 

The  Bachman-Landau  0  -  o  notation; 

Let  x  be  a  variable  element  of  a  topological  T^-space,  and  let 

it  range  over  a  set  R  ,  let  x^  be  a  cluster  point  of  R  (which  may  or  may 

not  belong  to  R) .  cp(x)  and  f(x)  denote  real-  or  complex- valued  functions 
defined  for  x  in  R  .  We  write 

(i)  cp  =  o(f)  in  R  ,  if  there  exists  a  constant  (i.e.  a  number 
independent  of  x)  A  such  that 

| Cp  |  <  A\f\  for  all  x  in  R  . 

(ii)  9  =  0(f)  as  x  x  ,  if  there  exists  a  constant  A  and  a 

neighbourhood  N. (x  )  of  x  such  that 
A* (i) (ii) * * v  o  o 

| Cp  |  <  A\lp\  for  all  x  in  RnN^(xQ) 

9  =  o(f)  as  x  ->  x^  ,  if  for  any  given  e>  0  there  exists  a 

N  (x  )  of  x  such  that 
e '  o  o 

1 9 1  <  e\f\  for  all  x  in  R^N£(x  ) 


(iii) 

neighbourhood 
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Consider  first 
N 

(»>  I  / 


V  P  (sv)  V 

o  n 


n-  o 


A+n- 1 


Ka-r)1 


(D5  1  !. 


+  V 


exp[-sv]  dv 


N  n 

I  I 

n=o  k=o 


(n)  1c 
Ck  S 


,v  A+n+k-1 


o  v 


J  o 
[D] 


(g-o* 

l-r2 


+  v 


exp[-sv]  dv 


By  equation  {j>h) ,  we  have 


(52) 


,k  r  o 


j 


o 

[D] 


A+n+k-1 


v 


(g-r); 

l-r2 


+  v 


exp[-sv]  dv 


,  p  oo  A+n+k-1 
k  /  v 


-'o  f^a-r)2 

^  1  l-r2  +V 


exp[-sv]  dv 


and  by  equation  (28) 


(53) 


*oo  A+n+k-1 

v 


■  °  i  <a-r)£  +  J 

I.  l-r2  J 


exp[-sv]  dv 


=  r(A+n+k) 


A+n+k- 1 

l  TT  (\_L-n  «_L-  .  .  si(a-r}f 


l-r'1 


U  (A+n  t-k  ;  A+n+k  j  ~A~~ ) 

l-r2 


Thus , 
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N  n 


<»>  1  1 


n=o  k=o 


(n)  k 

V  8 


o  A+n+k- 1 

v 


J  o 
[D] 


exp[-sv]  dv 


N  n 


* 


l  i  fef} 


A+n+k-1 


n=o  k=o 


U (A+n+k  )  A+n+k  ;  9  (a.~r )  )  , 

1-r2 


Now  the  order  of  magnitude  of 


n 


^  cW  (A+n+k) 
k=o 


A+n+k- 1 

r^-  T  U  (A+n+k  J  A+n+k  ;  s-(a"r^  ) 

2  J  1-r2 


is  uniformly  (in  s  only) 


0  (  TS  >  ’  as  s 


since  (c.  f.  equation  (29))  , 


(55) 


U  (A+n+k  •  A+n+k  ;  ).—  ) 

1-r2 


.  I  1=1 


^  (a-r) 


A+n+k 


pj  A+n+k 

s 


■l  +  0  (~  7  1  \p 

.  -r 
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Now  consider  the  remainder  term 


(56)  “  J 


,v  A- 1 

o  v 


[D] 


/^+v|  N 

t.  1-r2  1 


V  exp[-8v]  dv 


A+N 


-'o  f  (a-r)2  \ 

to)  \pr  +  v) 


(sv>v)  exp[-sv]  dv  , 


where  | f ^ |  <  M  ,  for  some  fixed  M  >  0  uniformly  for  |v|  <  2 |v^ 


Thus, 


(57) 


J 


r 


v  A+N 

o  v 


o  f  (a-r)2  \ 

t”)  \  1-r2  J 


f^(sv,v)  exp[-sv]  dv 


<  M 


P 

J  o 

[D] 


A+N 


v 


I1! 


-r  )2 

+ 


v 


exp[-sv]  dv 


<  M 


A+N 

r( A+n+1)  {  i  U  (A+n+1;  A+n+1;  P—P--  ) 

L  l-r2  1-r2 


Now  by  (29) 


■  .v;  ■:  !  '  '  ;•  '  -..i 


o>- 


(58) 


-  no  - 


U(A+N+1;  A+N+l; 

1-r2 


<  m 


Jihll 1 

l(a-r)2 


A+N+l 


1 

\  +  N  +  1 

s 


for  some  fixed  m  >  o  . 


Thus , 


(59) 


rv 

I  ° 


J 


o 

[D] 


a+n 


J  (a-r)2  1  N 


f..(sv,v)  exp[  -sv ]  dv 


<  mM 


t(A+n+i) 


(1-r2) 

(a-r)2 


X7TTTT 


Now  consider 


(60)  mM 


T (A+N+l) 


1  -  rg 
(a-r)2 


=  mM 


T (A+N+l) 


1  -  rg 
(a-r)2 


which  -*  0  as  s  ->  , 


Thus 


\  °  l  ~~X+N 


where  the  order  relation  is  uniform  in  s  but  not  in  r  , 
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Thus , 


(61) 


->  o 

[D] 


A-l 


fe?  *4 


G(sv,v)  exp[-sv]dv 


N  n 


I  I  -f>  «>»«  ft?} 


A+n+k-1 


n=o  k=o 


X  sk  U  ^A+n+lc  ;  A+n+k  ;  -s-^a~r  ^  ■■  ^ 
\  1-r2  ' 


=  o 


A+N 


which,  together  with  equation  (37)#  implies 


(62)  I*  = 


A-l  x 

V  f(v) 


J°  +  v} 

L  i_r2  ' 


exp[-sv  ^(v) ]  dv 


£  A  ,  s  r(q-r)2~l 

)  Ck"  r^+n+k)  l  1-r2  4 


A+n+k- 1 


n=o  k=o 


X  sku(Jv+n+k  ;  X+n+k  ;  =  o 


which  implies  that 


■  >  ,  i?  ;■  /> 


1 


exp[ -sv  lp(v)  ]  civ 


(6?) 


J  o 
[D] 


A-l  W  \ 

v  f(v) 


f  (a-r  Y 
\  1-r2 


+  v 


oo  n 


n=o  k=o 


^  f (A+n+k) 


A+n+k- 1 

f (a-r)2l 
L  l_r2  ^ 


k 

s 


/ 


A+n+k  •  A+n+k  ; 


s (g-r )2\ 
1-r2  / 


with  respect  to  the  sequence 


)  • 


Thus  from  (62)  we  have 


exp[ -sv  ^(v) ]  dv 


-  I 

n=r 


n 


> 


(n) 


T (A+n+k ) 


g-r ); 

l-r‘ 


2]  A+n+k  - 1 


k 

s 


U^A+n+k 


A+n+k 


s(a-r)^ 
1-r2  ' 


X 
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Thus, 

(64) 


J  o 
[D] 


A-l  ,{  s 
v f(v). 


g-r ); 

1-r2 


+  v 


exp[ -sv  ^(v) ]  dv 


N-l  n 


n=o  k=o 


C<n)  r(A+n+k){-^2 

1-r2 


A+n+k-1 


X  s^U/^A+n+k  ;  A+n+k  ;  S  ^~r ), 
\  1-r2 


N 

I 

k=o 


c  -N')  r(x+u+k)  J-(s~rl£ 

k  ll-r£ 


A+N+k-1 


sk  U  ^A+N+k;  A+N+kj  —  -- 


1- 


+  o 


A+N 


Now  by  (58), 


N 


(65)  ^  cf)  r(A+N+k) 

k=o  L 


A+N+k-1 

■r^1  sk  U/^A+N+k;  A+N+kj 


¥ 


l-r‘ 
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Thus, 

(66) 


A-l  \ 
v  f(v) 


exp[  -sv  ip(v)]  dv 


N 


n 


I  I 

n=o  k=o 


(n) 


T (A+n+k) 


W]£\ 

l-r2  1 


A+n+k- 1 


k 

s 


U(  A+n+k;  A+n+k; 


s(g-rn 

l-r2  J 


where  the  order  relation  is  uniform  in  s  but  not  in  r  . 


A  more  exact  result  for  the  estimate  of  the  order  of  magnitude 
of  the  error  involved  in  stopping  at  n  =  N  in  (63)  may  be  obtained  by 
use  of  equation  (k9)  .  Thus 


(67) 


iV 


J  o 

[D] 


A-l 


v 


(q-r)‘ 

l-r2 


+  v 


exp[-sv]  dv 


A-l 


J  o 
[D] 


f(a-r); 

'  l-rJ 


+  v 


0(vN+'*‘  exp[K|v 


sv 


]) 


exp[~sv]  dv 


[D] 


exp[ -sv{1-Kvq 


where  v  >  o  and  sv  >  o  for  all  v  on  D  in  [o,v  1  . 
o  L  7  o 
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Now  for  n  sufficiently  small  choice  of  v  ,  1  -  Kv  >  5>  0 

o  *  o  — 

for  some  fixed  6  .  Thus  for  some  fixed  €  >  o  ,  for  o  <  v  <  e  .  we 

o  — 


have 


and  since  8  ,  \  and  N  are  fixed  , 


=  0 


r(X+N+l)  U  (X+N+l  ;  A+N+l 


s(a-r)2  ^ 
1-r2 


(by  equation  (2^)). 
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Thus  we  have  shown 


(71) 


A-l  ,,  \ 
v  £(v) 


J  o 

[D] 


ISzll 

1-r2 


+  v 


exp[  -sv  V'(v)  ]  dv 


N 


n 


n=o  k=o 


A 

4n)  r(A+n+k)  {iSElfl 

k  l  1-r2  > 


A+n+k- 1 


k  . 
s  U 


^A+n+k; 


A+n+k; 


s  (g-r )2  \ 

1-r2  / 


+  0 


r  (A+n+1)  U (A+n+1 ; A+n+1  ; 


s(ci-r)2  ) 

1-r2 


Now  by  equations  ( 2 ),  (12)  and  (17)  we  have 


(72)  I (r ) 


- *.(£) -  (i.ts)s  at 

a-r-is/l-r^  t 


[D] 


v"2f1(v) 


(l-v)S  dv 


v2  f2(v) 


(l-v)S  dv 


where. 
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and. 


(73) 


and 


(7^) 


(i) 

V  =  t2 

(ii) 

fx(v) 

(iii) 

f2(v) 

2bfu? 


Cx(t)  -  x(-t)] 


(iv)  D  is  the  v-image  of  L*  . 
Further, 

V~2  f,  (v) 

(l-v)  dv  = 


f1 

v"2  fi(V) 

L 

[D] 

1  v  2  f^(v) 


J  o 
[D] 


v:>  f2(v) 


'  o 

CD] 


J 

l  1-r2 


+  v 


(l-v)S  dv  = 


J  o 
[D] 


Thus,  by  equations  (71)  and  (73) 
v'2  f1(v) 


N  n 


q-r  )* 

1-r2 


+  v 


exp[-sv  ( v ) ]  dv 


v2  fg(v) 


+  V 


exp[ -sv  f(v) ]  dv, 


(75) 


f  1 — -  (i-v)sdv=y  y 

Jo  f(a-r)2  1  L  L  k  1  1-r® 

[t>]  17^7- +  vj  n=°  k=0 

A  Sk  U  [^n+k+ljn+k+l; 

V  -I  „2 


+  0 


N+i 

r(N+  |)  u(H+  |  ;  n+  |  ;  .sfe-S)!)}  t 

1-r2  1-r2  / 
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and,  by  equations  (71)  and  (7k) 


(76) 


J  o 

[D] 


(g-r); 

1-r2 


+  v 


(l-v)S  dv 


N-l  n 


■  I  Z  -i"-81  «"*“  i'  fr?} 


n=o  k=o 


n+k+| 


k 

s  U 


f  ,  3  .  ,  3  .  s(oc-r)2  \ 

(n+k+  2  i  n+k+  2  ’  J 


N+| 

0  /T(SL-r.).2.|.  r(N+  |)  u(N+  |  ;  N+  |  ;  ^-a"r^ 
1-r2  J  1-r2 


where 

(n,l)  A 

c,  7  and 

k 

(n,2)  .  , 

c£  7  are  given  by 

(77) 

00  n 

I  z 

cO^1)  sk  vn+k  _  fi(v)  exp[ sv  (1  -  V'(v))  ] 

n=o  k=o 

and 

(78) 

00  n 

z  z 

c(n>2)  sk  vn+k  _  f g ( v )  exp[sv  [1  -  V'(v)]  ] 

n=o  k=o 

Now  changing  the  index  of  summation  in  (76)  and  (78)  we  have, 
corresponding  to  (76) 


-A’ 
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(79) 


1  v?  f  (v) 


K  1  (g-r)2  \ 

w  IT^~  +  v> 


(l-v)  dv 


N  n-1 


n=l  k=o 


n+k-i 


)  Z  r(n+k+5)  s''  U  ^n+k4;n+k4; 


+  o7J^=e1 

1-r2 


21  N+^ 


r(N+  |)  U(N+  |  ;  N+  |  ; 


where  c^n  is  given  by  (the  correspondent  to  ( 78 )) 


(80) 


00  n-l 

I  z 

n=l  k=o 


(n-l, 2)  k  n+k-1 
c.v  ’  8  v 

k 


f2(v)  exp[sv  [1  -  V'(v))  ] 


Thus  (72)  becomes 


(81) 


r: 

[L*] 


JLdI 


a-r-Wl-r2  t 


(l-t2)S  dt 


N  n 


■  S.  Z  Z 


lY 


n+k~|- 


n=o  k=o 


s  v  U  (  n+k+i-;  n+k+4, 

\  '  1-r2 


N  n-l 


-  Z  Z 


2ln+k'iSk  D  In+k+fen+k+i- 


n=l  k=o 


01  v  -  -  2>  l,TRT2>  o 

2J  k  l~r2 


N+| 


+  °({^}  r(H+|)U(N+|;N+|;^£)) 


APPENDIX  IV 


GRAPHICAL  RESULTS 


In  this  appendix  seven  graphs  are  given.  The  purpose  of 
drawing  these  graphs  is  to  indicate  the  general  trend  and  shape  of 
the  distributions  obtained  in  Chapter  three:  specifically,  the  results 
given  in  equations  (3.2.14)  and  (3.2. 15)  are  for  the  circular  case,  and 
in  equations  (3.2. 16)  and  (3.2. I7)  for  the  stationary  case. 

The  seven  graphs  are  divided  into  three  groups.  Each  group 
is  preceded  by  a  set  of  prefatory  remarks  regarding  the  characteristic 
features  brought  into  prominence  by  that  particular  group. 

As  regards  the  graphs  themselves,  the  variable  r  varies  over 
the  closed  interval  [-l,+l],  and  is  represented  along  the  horizontal 
axis.  In  Groups  A  and  B,  h(r)  and  H(r)  [Daniels'  density  function 
and  our  distribution  function  respectively]  are  represented  along  the 
vertical  axis.  The  scale  for  h(r)  being  given  on  the  left  and  for 
H(r)  being  given  on  the  right.  In  group  C  only  the  distribution 
function  H(r)  is  represented  along  the  vertical  axis. 
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GROUP  A. 


Each  of  the  following  three  graphs  contains  a  tracing  of  the 
function  given  by  equations  (3.2.14)  and  (3.2.15),  which  is  the  approximate 
Distribution  function  for  the  circular  sample  serial  correlation  coefficient. 
Included  also  is  a  tracing  of  the  function 

n  1 

n 

(l-2or+a2)2 

which  is  the  approximate  Density  function  for  the  same  sample  serial 
correlation  coefficient  as  derived  in  DANIELS  (1956). 

The  three  cases  considered  are 


(i) 

n  =  10, 

a  = 

-0.5  . 

(ii) 

n  =  20, 

a  = 

0.0  . 

(iii) 

n  =  30, 

a  = 

+0 .5  • 

These  three  cases  give  a  fair  cross-section  of  the  behavior  of  the 
derived  approximate  distribution  function,  and  its  relationship  with 
Daniels'  result. 

The  most  prominent  feature  exhibited  by  the  three  graphs  is 
that  the  skewness  of  the  distribution  in  question  is  a  function  of  a  . 


.rJjI+l) 


r(§  +  |) 
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Distribution  and  Density  Functions  for  n  =  20,  a-  =  0.0 
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GROUP  B . 


The  following  two  graphs,  also  for  the  circular  case,  show  the 
effect  on  the  Kurtosis  of  distribution  in  question  by  changes  in  the 
size  of  n.  The  two  cases  considered  are 

(i)  a  =  0.0,  n  =  10 

(ii)  a  =  0.0,  n  =  20  . 

The  graphs  indicate  that  the  Kurtosis  is  a  function  of  n. 

Both  group  A  and  group  B  serve  to  indicate  the  relationship 
between  Daniels'  result  and  ours.  Although  the  relationship  is  approx¬ 
imate,  there  seems  to  be  an  excessive  discrepancy  near  r  =  -1,  This 
discrepancy  we  attribute  to  the  fact  that  neither  Daniels'  result  nor 
ours  is  ~eliable  in  the  regions  near  r  =  +  1. 

A  similar  relationship  between  the  Distribution  function  and 
Daniels'  Density  function  exists  for  the  stationary  case.  Because  of 
this  similarity  between  the  stationary  and  circular  cases  we  have 
omitted  the  sketches  (as  in  groups  A  and  B)  for  the  stationary  case. 
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GROUP  C. 


For  the  stationary  case,  the  following  graphs  indicate 

(i)  for  fixed  n,  the  variation  in  the  distribution  function 
for  three  different  values  of  a. 

(ii)  for  fixed  a,  the  variation  in  the  distribution  function 
for  three  different  values  of  n. 

To  bring  into  focus  the  obvious  variations  a  system  of  super¬ 
position  has  been  adopted,  the  various  cases  considered  are  clearly 
marked  on  the  respective  graphs. 

The  results  of  these  two  graphs  overlap  somewhat  with  those 


given  by  groups  A  and  B. 
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Fig.  C. 1  Distribution  Functions  for 


i.  ■' 


100 
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Fig.  C.2  Distribution  Functions  for  a  =  0.0 


